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1 TeV Astronomy 

1.1 Introduction 

The discovery of more than 100 extraterrestrial sources 
of Very High Energy (VHE, ^ 100 GeV) or TeV gamma- 
radiation belongs to the most remarkable achievements 
of the last decade in astrophysics. The strong impact 



of these discoveries on several topical areas of modern 
astrophysics and cosmology are recognised and highly 
appreciated by different astronomical communities. The 
implications of the results obtained with ground-based 
TeV gamma-ray detectors are vast; they extend from 
the origin of cosmic rays to the origin of Dark Mat- 
ter, from processes of acceleration of particles by strong 
shock waves to the magnetohydrodynamics of relativistic 
outflows, from distribution of atomic and molecular gas 
in the Interstellar Medium to the intergalactic radiation 
and magnetic fields. 

TeV gamma-rays are copiously produced in environ- 
ments where effective acceleration of particles (electrons, 
protons, and nuclei) is accompanied by their inten- 
sive interactions with the surrounding gas and radiation 
fields. These interactions contribute significantly to the 
bolometric luminosities of young Supernova Remnants 
(SNRs), Star Forming Regions (SFRs), Giant Molecular 
Clouds (GMCs), Pulsar Wind Nebulae (PWNe), com- 
pact Binary Systems, Active Galactic Nuclei (AGNs) and 
Radio Galaxies (RGs), etc.. 

The fast emergence of gamma-ray astronomy from an 
underdeveloped branch of cosmic-ray studies to a truly 
astronomical discipline is explained by the successful re- 
alization of the great potential of stereoscopic arrays 
of Imaging Atmospheric Cherenkov Telescopes (IACTs) 
which act as effective multifunctional tools for deep stud- 
ies of cosmic gamma-radiation. 

Being recognised as one of the most informative win- 
dows to the non- thermal Universe, the VHE domain of 
gamma-rays provides also means for probing fundamen- 
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tal physics beyond the reach of terrestrial accelerators. 
In particular, the indirect search for Dark Matter and 
tests of quantum gravity using these energetic gamma- 
rays are amongst the high priority objectives of the cur- 
rent and future projects with involvement of ground- 
based gamma-ray detectors. In this regard, the TeV 
gamma-ray astronomy is considered a key component 
of the new interdisciplinary research area called Astro- 
Particle Physics. 

1.2 Ground-based Detection Technique 

The atmosphere of the Earth is not transparent to 
gamma-rays, therefore their direct registration requires 
platforms in space. The currently operating Fermi Large 
Area Telescope (Fermi-L AT; formerly GLAST) is a pow- 
erful satellite-borne instrument designed for deep sur- 
veys with a very large field view of the order of 2 stera- 
dian. Presently, the study of the sky in MeV and GeV 
gamma-rays by Fermi-LAT is complemented by a some- 
what smaller-scale telescope on the italian X-ray and 
gamma-ray satellite AGILE (Astro-rivelatore Gamma a 
Immagini LEggero). The angular resolution of these in- 
struments below 1 GeV is quite modest (larger than 1°), 
but it becomes significantly better at higher energies, ap- 
proaching 0.1° above 10 GeV. Fermi-LAT covers a very 
broad energy region of primary gamma-rays extending 
from tens of MeV to hundreds of GeV (HE; up to 300 
GeV). However, beyond 10 GeV the gamma-ray fluxes 
are generally very faint, so that the effective detection 
area of Fermi-LAT cannot provide adequate statistics 
for comprehensive spectral and temporal studies in the 
VHE domain. 

There is not much hope that space platforms could of- 
fer, in any time in the foreseeable future, instruments 
with detection areas significantly exceeding 1 m 2 . This 
dramatically reduces the potential of studies of VHE 
gamma-rays from space. Fortunately, at these energies 
an alternative method can be exploited based on the 
registration of atmospheric showers, either directly or 
through their Cherenkov radiation. The faint and brief 
Cherenkov signal of relativistic electrons produced dur- 
ing the development of the electromagnetic cascades in 
the atmosphere, lasts only several nanoseconds, but it 
is sufficient for detection using large optical reflectors 
equipped with fast optical receivers. With a telescope 
consisting of a 10 m diameter reflector and a multichan- 
nel camera of pixel size of ~ 1/4 degree and a field- 
of-view of ~ 3 degree, primary gamma-rays of energy 
E ^ 100 GeV can be effectively collected across ground- 
level distances as large as 100 m providing a huge area 
for the detection of primary gamma-rays, A e R ^ 10 4 m 2 . 
The total number of photons in the registered Cherenkov 
light image is proportional to the primary (absorbed in 
the atmosphere) energy, the orientation of the image cor- 



relates with the arrival direction of the gamma-ray pho- 
ton, and the shape of the image contains information 
about the origin of the primary particle (a proton or 
a photon?). The stereoscopic observations of air show- 
ers with two or more telescopes located at distances of 
about 100 m from each other, provide effective rejection 
of hadronic showers (by a factor of 100), as well as good 
angular resolution (better than 0.1°) and energy reso- 
lution (better than 15 per cent). At energies around 1 
TeV, this results in a minimum detectable energy flux as 
low as 3 x 10 -13 erg/cm 2 s (see e.g. [1]) This is much bet- 
ter than in any other gamma-ray domain, including the 
GeV energy band, where the sensitivity of Fermi LAT 
cannot compete with the performance already achieved 
by the H.E.S.S., MAGIC and VERITAS IACT arrays 
in the TeV energy band. Thanks to the very large col- 
lection area, the IACT technique provides large gamma- 
ray photon statistics even from relatively modest TeV 
gamma-ray emitters. In combination with good energy 
and angular resolutions, the gamma-ray photon statistics 
appears to be adequate for deep morphological, spectro- 
scopic and temporal studies. This also makes the IACT 
arrays powerful multifunctional and multi-purpose tools 
for the exploration of a broad range of non-thermal ob- 
jects and phenomena. The potential of the IACT arrays 
has been convincingly demonstrated by the H.E.S.S., 
MAGIC and VERITAS collaborations (see, e.g. [2], and 
references therein). 

The IACT arrays are capable to study not only point- 
like, but also extended sources with an angular size up 
to 1 degree or somewhat larger. Moreover, the high 
flux sensitivity and relatively large (^4°) field of view 
of IACT arrays allow rather effective all-sky surveys as 
demonstrated by the H.E.S.S. collaboration. On the 
other hand, the potential of IACT arrays is rather lim- 
ited for the search of very extended structures like the 
galactic plane diffuse emission or the huge radio lobes of 
the nearby radio galaxy Centaurus A. IACT arrays have 
a limited capability for "hunting" of solitary events like 
the possible VHE counterparts of Gamma Ray Bursts. 
In this regard, the detection technique based on direct 
registration of particles (leptons, hadrons and photons) 
of extensive air showers (EAS) is a complementary ap- 
proach to the IACT technique. 

The traditional EAS technique, based on scintillators or 
water Cherenkov detectors spread over large areas, was 
originally designed for the detection of cosmic rays at 
PeV and EeV energies. In order to adopt this tech- 
nique to gamma-ray astronomy, the energy threshold 
needs to be reduced by two or three orders of magni- 
tude. This can be achieved using dense particle arrays 
located on very high altitudes. The feasibility of both 
approaches recently have been successfully demonstrated 
by the ARGO and Milagro collaborations. In particular, 
several very extended sources have been reported by the 
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Milagro group. These results, as well as the potential for 
continuous monitoring of a significant part of the sky, 
which might lead to exciting discoveries of yet unknown 
VHE transient phenomena in the Universe, strongly sup- 
port the proposals of constructing high altitude EAS de- 
tectors (see [1] for a review) like HAWC, a High Altitude 
Water Cherenkov Experiment, presently under construc- 
tion on a site close to Sierra Negra, Mexico [3]. The 
5 yr-survey sensitivity of HAWC above 1 TeV is expected 
to be comparable to the sensitivity of Fermi-LAT at 1 
GeV. Thus HAWC will be complementary to Fermi for 
continuous monitoring of more than 1 steradian fraction 
of the sky at TeV energies. At higher energies, recently 
a new project called LHAASO (Large High Altitude Air 
Shower Observatory) has been suggested. The proposed 
huge detector facility at Yangbajing, China, will consist 
of several sub-arrays for the detection of the electromag- 
netic and muon components of air showers. They will 
cover a huge area, and can achieve an impressive sen- 
sitivity at energies of several tens of TeV (see Fig. 1). 




^ Q-15 i | i | i I i I i I 
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E [GeV] 

Fig. 1 Energy-flux sensitivities of current and future ground- 
based detectors - the IACT and EAS arrays in the energy range 
10 10 to 10 15 eV (courtesy of G. Sinnis). 

1.3 Future IACT Arrays 

The future of observational gamma-ray astronomy, at 
least for the next 10-15 years, is connected with the next- 
generation IACT arrays, first of all with the observatory 
CTA (Cherenkov Telescope Array) [4], cf. also Fig. 2. 
The next generation of IACT arrays are aiming at (i) a 
significant (by an order of magnitude) improvement of 
the flux sensitivities in the standard 0.1-10 TeV energy 
interval ( Te V regime) , and (ii) an expansion of the en- 
ergy domain of IACT arrays in both directions - down 
to 10 GeV (multi-GeV regime) and well beyond 10 TeV 
(sub-PeV regime): 




Fig. 2 Possible layout of the next-generation CTA instrument. 
From Ref. [4]. 

• TeV regime: 

This is the "nominal" energy region where the IACT 
technique has achieved its best performance. The po- 
tential in this energy regime is still not saturated. With 
a stereoscopic array consisting of tens of 10 m-diameter 
(medium-size) class telescopes the minimum detectable 
energy flux could be reduced to the level of 10 -14 
erg/cm 2 s, and the angular resolution be improved to 
66 ^ 3 arc minutes. Generally, the optimum distance 
between the telescopes is considered to be around 100 m, 
the radius within which the Cherenkov light is dis- 
tributed more or less homogeneously. However, if highest 
priority is given to the performance at energies around 
1 TeV and beyond, an increase of the distance between 
telescopes up to 300 m could be an attractive option. 
For a fixed number of telescopes this would increase 
the detection area by an order of magnitude, and, at 
the same time, improve the angular resolution to 1-2 
arc minutes, although at the expense of a somewhat 
higher (by a factor of two or three) energy threshold. 
In any case, a reduction of the minimum detectable en- 
ergy flux around 1 TeV down to 10 -14 erg/cm 2 s seems 
to be a challenging but feasible "target". It will be a 
great achievement even by the standards of the most ad- 
vanced branches of observational astronomy, allowing us 
to probe, in particular, potential TeV gamma-ray sources 
at luminosity levels of 10 32 (d/10 kpc) 2 erg/s for galac- 
tic and 10 40 (d/100 Mpc) 2 erg/s for extragalactic ob- 
jects. Although for moderately extended sources, e.g. 
of angular size ^ ~ 1°, the minimum detectable energy 
flux will be by a factor of ~ 10 — 30 higher, it 

would compete or be better than the energy flux sen- 
sitivities of the best current X-ray satellites, Chandra, 
XMM-Newton, INTEGRAL and Suzaku, and thus allow 
the deepest probes of non-thermal high energy phenom- 
ena in extended sources, in particular in shell- type SNRs, 
Giant Molecular Clouds, Pulsar-driven Nebulae (Pleri- 
ons), Clusters of Galaxies, hypothetical Giant Pair Halos 
around AGN, etc. Such a system of 10-12 m diameter 
class IACTs with a field of view (FoV) of 6-8 degrees, 
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will most likely constitute the core of the Cherenkov 
Telescope Array (CTA) - an initiative towards a major 
ground-based gamma-ray detector (see Fig. 2). 

• Sub-PeV regime: 

External and intergalactic absorption of gamma-rays, 
the limited efficiency of particle acceleration, the escape 
of highest energy particles from the source etc., result in 
a suppression of fluxes at the highest energies. The gen- 
eral tendency of decreasing gamma-ray fluxes with en- 
ergy becomes especially dramatic above 10 TeV. There- 
fore, any meaningful study of cosmic gamma-rays be- 
yond 10 TeV typically requires detection areas as large 
as 1 km 2 . An effective and straightforward approach 
would be the use of IACT arrays optimised for detection 
of gamma-rays in the region up to 100 TeV and beyond. 
This can be realised by modest, approximately 10-30 m 2 - 
area reflectors separated from each other, depending on 
the scientific objectives and the configuration of the im- 
agers, by 300 to 500 m. The requirement on the pixel 
size of imagers is also rather modest, 0.25° or so, how- 
ever they should have large, up to 10 degree FoV for 
simultaneous detection of showers from distances more 
than 300 m [5]. A sub-array consisting of several tens of 
such small-size telescopes is included in the concept of 
CTA with a primary goal to study the energy spectra of 
gamma-ray sources well beyond 10 TeV. It will serve as 
a powerful tool for searches of galactic cosmic ray "Pe- 
Vatrons", as well as nearby (R <C 10 Mpc) radio and 
star burst galaxies. 

• Sub-100 GeV regime: 

The energy threshold of detectors, £^ t h, is generally de- 
fined as a characteristic energy at which the gamma-ray 
detection rate for a primary power-law spectrum with a 
photon index 2-3 achieves its maximum. It is known from 
Monte Carlo simulations as well as from the experience of 
operation of previous generation of IACTs, that in prac- 
tice the best sensitivity is achieved at energies exceeding 
several times E^. Thus, for optimisation of gamma-ray 
detection around 100 GeV, one should reduce the en- 
ergy threshold of telescopes to E t h ^ 30 GeV. This can 
be done by using very large, 20 m-diameter (large-size) 
class reflectors. On the other hand, the reduction of the 
threshold to 30 GeV is an important scientific issue in 
its own right; the intermediate interval between 30 and 
300 GeV could be crucial for certain classes of galactic 
and extragalactic gamma-ray sources. A sub- array con- 
sisting of several large-size telescopes as foreseen in CTA 
(see Fig. 2) will indeed significantly broaden the topics 
and scientific objectives of CTA. 

Each of the IACT arrays discussed above covers at least 
two decades in energy with significant overlaps of the 
energy domains. Since these arrays contain the same 



basic elements, and generally have also common scien- 
tific motivations, an ideal arrangement would be if these 
sub-arrays are combined in a single facility which would 
have a sensitive and homogeneous coverage throughout 
the energy region from approximately 30 GeV to 300 
TeV. The concept of CTA is based, to a large extent, 
on this argument [4]. The high detection rates, coupled 
with good angular and energy resolutions over four en- 
ergy decades will make CTA a powerful multi-purpose 
gamma-ray observatory with a great capability for spec- 
trometric, morphological and temporal studies of a di- 
verse range of persistent and transient high-energy phe- 
nomena in the Universe. 

• Mult i- GeV regime: Gamma-Ray Timing Explorers 

Despite the recent great achievements of high energy 
(HE) gamma-ray astronomy, there are obvious shortcom- 
ings in the performance of the current so-called "pair- 
conversion" tracking detection technique - the most ef- 
fective approach used in satellite-borne instruments for 
detection of gamma-rays at energies above 100 MeV. One 
should note that the flux sensitivity of Fermi-LAT at 1 
GeV of about 10 -12 erg/cm 2 s can be achieved only after 
one year all-sky survey. While for persistent gamma-ray 
sources this seems to be an adequate sensitivity (given 
that a huge number of sources are simultaneously moni- 
tored within the large and homogeneous FoV), the small 
detection area significantly limits its potential, in par- 
ticular for detailed studies of the temporal and spectral 
characteristics of highly variable sources like blazars or 
solitary events like gamma-ray bursts (GRBs). It will not 
be easy to improve the sensitivity achieved by Fermi- 
LAT at high energies by any future space-based mis- 
sion, unless the Moon would be used in the (far) fu- 
ture as a possible platform for an installation of very 
large (^> 10m 2 ) area pair-conversion tracking detectors. 
Apparently, the space-based resources of GeV gamma- 
ray astronomy have achieved a point where any further 
progress would appear extremely difficult and very ex- 
pensive. In any case, for the next decades to come there 
is no space-based mission planned for the exploration of 
the high-energy gamma-ray sky. On the other hand, the 
principal possibility of an extension of the IACT tech- 
nique towards 10 GeV promises a new breakthrough in 
gamma-ray astronomy [1]. The (relatively) large gamma- 
ray fluxes in this energy interval, together with the huge 
detection areas offered by the IACT technique, can pro- 
vide the highest gamma-ray photon statistics compared 
to any other energy band of cosmic gamma-radiation. 
Thus, in the case of a realization of 10 GeV-threshold 
IACT arrays, the presently poorly explored interval be- 
tween 10 and 100 GeV could become one of the most 
advanced domains of gamma-ray astronomy with a great 
potential for the studies of highly variable phenomena. 
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The reduction of the energy threshold down to 10 GeV 
or even less is principally possible within the basic con- 
cept of the I ACT technique, but it requires an extreme 
approach of using 25 m diameter class telescopes with 
very high 40%) quantum efficiency focal plane im- 
agers to operate in a robotic regime at very high (5 km 
or) mountain altitudes [6]. 

The energy range from several GeV to 30 GeV has very 
specific astrophysical and cosmological objective: explo- 
ration of the highly variable non-thermal phenomena in 
the remote universe at redshifts of z = 5 (like large red- 
shift quasars and GRBs) , as well as the study of compact 
galactic sources such as pulsars and microquasars. A re- 
alization of such a gamma-ray timing explorer, hopefully 
during the lifetime of the Fermi observatory would be a 
great achievement for gamma-ray astronomy. 



2 TeV Sources 

2.1 Supernova Remnants 

Massive stars are believed to end their life undergoing a 
supernova explosion. This explosion blows off their other 
layers into a supernova remnant (SNR), which heats the 
surrounding medium and accelerates cosmic-rays (elec- 
trons and protons) to extremely high energies. The ra- 
diation from shell-like SNRs consists of thermal emis- 
sion from shock-heated gas and non-thermal emission 
from shock-accelerated particles. The theory of diffusive 
shock acceleration (DSA) at shock fronts [7,8] predicts 
the production of a population of accelerated particles 
in SNRs that can interact with ambient magnetic fields, 
with ambient photon fields, or with matter. The amount 
of relativistic particles increases with time as the SNR 
passes through its free expansion phase, and reaches a 
maximum in the early stages of the Sedov phase. Corre- 
spondingly, the peak in gamma-ray luminosity typically 
appears some 10 3 -10 4 years after the supernova explo- 
sion. 

In the TeV domain, presently seven shell-type SNRs - 
Cas A [9-11], Tycho [12], SN 1006 [13], RX J1713.7-3946 
[14,15], RX J0852-4622 (Vela Junior) [16], RCW 86 [17], 
and G353.6-0.7 (HESS J1731-347)[18] have been firmly 
identified as VHE gamma-ray emitters (see Table 1). 
Remarkably, while the first six sources are well estab- 
lished young SNRs, the object G353.6-0.7 is the first 
SNR discovered serendipitously in TeV gamma-rays, 
and only later confirmed by radio and X-ray observa- 
tions [19,20]. Moreover, a possible new SNR candi- 
date, HESS J1912+101, has been postulated recently [21] 
based solely on its shell-type morphology at TeV ener- 
gies, although no counterpart at lower energies has been 
detected so far. The two latest examples demonstrate 



the potential of large field-of-view Cherenkov telescopes 
for serendipitously discovering extended SNRs (of typical 
size 0.2-1° at a distance up to ~3.5 kpc). Their relatively 
large sizes and 7-ray luminosities of about (0.1 — 1) x 10 33 
erg/s have enabled the detection of these objects up to 
distances of ~ 3.5 kpc (e.g., Tycho) with current instru- 
ment sensitivities (cf. Fig. 1). If the VHE gamma-ray 
luminosities detected from these objects reflect the typi- 
cal luminosity of the SNR population in the Galaxy, fu- 
ture instrument like CTA should be able to detect SNRs 
up to 15 kpc, thus sampling the whole Galaxy. Tak- 
ing the spatial distribution of SNRs in the Galaxy, their 
explosion rate, and the duration of the TeV emission (be- 
lieved to last a few thousand years) into account, roughly 
^100 new SNRs could be discovered at TeV energies [22] 
(in a naive approximation, without considering energy 
cut-offs, hard/soft spectral indices, etc.). Such an en- 
larged population would allow the study of these objects 
at different evolutionary stages, sampling their spectral 
energy distribution from a few hundred of MeV (with 
Fermi-LAT and AGILE) up to 100 TeV, in the cut-off 
regime. 
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**> 
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0-2 dee 


SN1006 


G353.6-0.7 
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1X2 dag 







Fig. 3 Example of four shell-type SNRs detected at TeV energies 
with the H.E.S.S. instrument. 

The sizes of several of these shell-like SNRs (> 0.1°) 
has allowed to resolve them in VHE (see Fig. 3). The 
images of SNRs such as SN 1006, RXJ1713. 7-3946 or 
RXJ0852-4622 have revealed a good correlation of the 
TeV emission sites with the non-thermal emission de- 
tected in X-rays, probing acceleration of relativistic par- 
ticles up to mufti- TeV energies. However, the relative 
contributions of accelerated protons and electrons to the 
gamma-ray production still remain unknown. The prob- 
lem is that the ratio of gamma-rays produced by ac- 
celerated protons interacting with the surrounding gas, 
and by ultra-relativistic electrons up-scattering the 2.7K 
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Table 1: Shell-like SNRs firmly detected at TeV energies 



Name 


Dist (kpc) 


Size (pc) 


Age (yrs) 


L 7 (10 33 erg/s) 


r 


RXJ1713.7-3946 


1 


17.4 


1.6 


8 


2.0 


RXJ0852-4622 


0.2(1) 


6.8(34) 


0.4(5) 


0.26(6.4) 


2.2 


RCW86 


1(2.5) 


11(28) 


1.6(10) 


1(6) 


2.5 


SN 1006 


2.2 


18.3 


1 


1.24 


2.3 


Cas A 


3.4 


2.5 


350 


7 


2.4 


Tycho 


3.5 


6 


438 


0.1 


1.95 


SNR G353.6-0.7 


3.2 


27 


2.5(14) 


10 


2.3 



CMB radiation, is very sensitive to generally unknown 
parameters, in particular to the gas density and the 
magnetic field of the ambient medium (cf., e.g. [302]). 
The efficiency of inverse Compton (IC) scattering is es- 
pecially high at TeV energies (up to E e ^100 TeV, it 
proceeds in the Thomson regime, with a corresponding 
cooling time t^ ol oc 1/E e oc 1/E^ 2 ). For example, the 
typical production time of a 1 TeV-photon by an elec- 
tron and a proton of the same characteristic energy of 
about 20 TeV, are ^5xl0 4 yr and 5xl0 7 (n/l cm 3 )- 1 
yr, respectively (see, e.g. [23]). Correspondingly, at 
1 TeV the ratio of the production rates of IC to 7r°- 
decay gamma-rays, is approximately 10 3 (W e /W p )(n/1 
cm -3 ) -1 , where W e and W p are the total energies in 
20 TeV electrons and protons, respectively. Thus, even 
for a very small electron-to-proton ratio (at the stage 
of acceleration), e/p = 10 -3 , the contribution of the IC 
component will dominate over the 7r°-decay gamma-rays 
(in the shell with a typical gas density n^lcm -3 ), unless 
the magnetic field in the shell significantly exceeds 10/iG. 
In this case, the accelerated electrons are cooled predom- 
inantly via synchrotron radiation, thus only a small frac- 
tion, wqmb/^b ~ 0.1(B/10/iG) -2 , will be released in IC 
gamma-rays. Alternatively, the proton-to-electron accel- 
eration ratio should exceed e/p ~ 10 3 which, in principle, 
cannot be excluded given the uncertainty associated with 
one of the key aspects of DSA related to the so-called in- 
jection problem (see [24]). 

In cases like RXJ1713. 7-3946, Tycho or Cas A, the mag- 
netic field has been estimated from multi- wavelength ob- 
servations to be >0.1 mG [25,26], restricting the contri- 
bution of the IC emission and in principle favouring an 
hadronic origin of the TeV emission. Nevertheless, if 
the IC and synchrotron components of the radiation are 
formed in different zones, these constraints are less ro- 
bust. For instance, a difference of the magnetic field in 
the upstream and the downstream region could result in 
a positional shift of the production regions of synchrotron 
X-rays and IC gamma-rays, and more complex models 
implying multi-zone emission would need to be invoked 
[27,28]. In general, while the distribution of the X-ray 



radiation is dominated by the strength of the magnetic 
field, the TeV emission traces the particle distribution 
and does not depend on the magnetic field, allowing a 
more unbiased study of the particle acceleration in the 
shell. With the angular resolution of current instruments 
(of the order of «0.1°) those different sites are still indis- 
tinct, but the future improvement of the angular resolu- 
tion to a few arcmin should permit a detailed study of 
the TeV radial profile in sources like RXJ1713. 7-3946 or 
SN 1006 in comparison with the X-ray radiation profile. 
The spectral energy distribution (SED) of these young 
SNRs extends over almost five decades, from a few hun- 
dred MeV to a few tens of TeV. At low energies the SED 
part for some of these TeV shell-like SNRs has been de- 
tected with the Fermi-L AT telescope [29-32]. The cov- 
erage of the spectrum at low energies has improved our 
understanding of the origin of the gamma-ray emission, 
but also evidenced a more complicated scenario in which 
different regions can contribute to the total emission, 
such as the reverse shock [28] or dense clouds embed- 
ded in the shock [33,34]. The photon spectra of Tycho, 
RXJ0852-4622 and Cas A continues to the MeV-GeV 
range with a rather hard spectral index of ~2.0 as pre- 
dicted by the DSA theory [35-37]. This fact, together 
with the high magnetic field amplification derived from 
synchrotron X-ray filaments, preventing in principle a 
large IC contribution from leptons, favour an hadronic 
scenario in these SNRs. Moreover, high-energy radia- 
tion up to at least a few TeV has also been detected 
from these SNRs without an indication for a turnover 
in the spectrum. An extension of the high-energy emis- 
sion by a factor two or three beyond 10 TeV could only 
be explained through hadronic interactions, given the 
fast Klein-Nishina-cooling suffered by 100 TeV-electrons 
emitting in this energy regime, and would robustly ex- 
clude an IC origin of the radiation. It would also provide 
a definitive probe of SNRs as origin of the cosmic-ray sea 
(see Section 3.1.). 

RXJ0852-4622 and RX J 171 3. 7-3946, for which large 
magnetic fields have been estimated, face some difficul- 
ties when modeling their gamma-ray emission. These 



Frank M. RIEGER; Emma de ON A- WILHELMI, and Felix A. AHARONIAN, Front. Phys. 



7 



two SNRs have similar ages, sizes, and radio, X-ray 
and TeV gamma-ray spectra, although RXJ1713. 7-3946 
shows a softer spectral index (~1.5) in the 100 MeV to 
1 GeV-energy range, similar to the predicted indices in 
a leptonic scenario. In both cases, the apparent low 
gas density (n~0.1 cm -3 ) [38] poses troubles to stan- 
dard hadronic scenarios [28,39,40]. Still, even in the case 
of a very low gas density of the shell, the contribution 
of hadronic gamma-rays could be significant, if acceler- 
ated protons interact with the dense cores of molecu- 
lar clouds embedded in the shell [82]. In such a case, 
slow diffusion could prevent low-energy particles to pen- 
etrate into these dense cores, suppressing the low-energy 
gamma-ray emission and naturally explaining the hard 
gamma-ray spectrum measured in RX J 1713. 7-3946. At 
the highest energies, RX J 171 3. 7-3946 shows a energy 
cut-off above few TeV, excluding PeV protons from this 
remnant. However, an escape of high-energy protons 
that cannot be confined in the shell, can not be excluded 
and might be a plausible explanation. In fact, at GeV 
energies, a large number of mid- age SNRs has been dis- 
covered, while only a small fraction of them shines at 
TeV energies. The gamma-ray emission in these cases is 
likely related to interactions of cosmic-rays with dense 
gaseous complexes [34]. In cases like W51C, detected 
up to ~5 TeV [41,42], an enhancement of the hadronic 
origin due to the large gas density in the region seems 
clearly favoured. On the other hand, the best example 
illustrating the escape of high-energy particles is the 10 4 
yr-old SNR W28 [43] , where a clear correlation between 
the TeV emission and massive molecular clouds emitting 
in CO has been observed. Some of these clouds are also 
bright at GeV energies. Another example of this type 
of scenario is IC 443 [44-47], where the GeV and TeV 
emission appear shifted from each other. These images 
seem to support an escape scenario where, depending on 
the location of the massive clouds, the time of particle 
injection into the interstellar medium and the diffusion 
coefficient, a broad variety of energy distributions may 
be expected. 

2.2 Pulsars 

Pulsars - rapidly rotating and highly magnetised neu- 
tron stars surrounded by a rotating magnetosphere and 
accompanied by relativistic outflows - emit radiation 
at all wavelengths. Charged particles (electrons and 
positrons) are thought to be efficiently accelerated in 
the electromagnetic fields of the pulsar, producing 7- 
radiation via e.g. curvature processes and supporting the 
formation of a cold relativistic outflow beyond the light 
cylinder. This pulsar wind carries almost the entire ro- 
tational energy of the pulsar in the form of Poynting flux 
and/or kinetic energy of the bulk motion, and creates a 
standing shock wave (the termination shock) when it in- 



teracts with the ambient medium. Particles accelerated 
at this shock are responsible for the steady and usually 
very extended non-thermal radiation observed (see Sec. 
2.3). 

Although pulsars have been traditionally a subject of 
radio astronomy, with ^1800 pulsars found beaming ra- 
dio waves, most of their radiation is emitted at high- 
energies (a few percent of their spin-down power). In- 
deed, in the last three years, the number of gamma-ray 
pulsars has increased exponentially from half a dozen to 
more than 150 [48] thanks to the new sensitive instru- 
ments Fermi-LAT and AGILE. Despite the high Galactic 
background, the periodic gamma-ray emission stands out 
due to the high fluxes, hard spectral index and power- 
ful timing analysis tools. The large statistics and good 
data quality has provided new insights into the physics 
of pulsars. In general, it is believed that the pulsed, 
periodic gamma-ray radiation originates in regions of 
the magnetosphere, called gaps, where the electric field 
has a parallel component along the magnetic field lines. 
This electric field efficiently accelerates electrons and 
positrons to relativistic energies causing them to emit 
synchro-curvature radiation in the form of gamma-rays. 
There are currently a few models that differ, primarily, 
on the location of these gaps [49-51], which are capable 
to explain the light-curves and spectral energy distribu- 
tions. Other mechanisms have also been suggested such 
as a magnetosphere with a force-free structure [52] or a 
striped wind topology [53]. The Fermi-LAT- measured 
light curves and energy spectra indicate that gamma-ray 
emission from the brightest pulsars is produced in the 
outer magnetosphere with fan-like beams scanning over 
a large portion of the celestial sphere. The energy spec- 
tra for most of the gamma-ray pulsars are best described 
by a power-law function with an exponential cutoff of 
the form E _r exp [— (E/E ) b ] with 6^1, and cut-off en- 
ergy E between 1 and 10 GeV [48]. The detection of 
gamma-rays beyond a few GeV without indication for 
a super-exponential attenuation (i.e., b > 1) effectively 
excludes the so-called polar cap model and gives a pref- 
erence to models of gamma-ray production in the outer 
magnetosphere (in order to avoid severe pair-production 
in the strong magnetic field in low-altitude zones). Most 
of the measured spectra can be well-fitted with a sim- 
ple exponential attenuation (6 = 1) [48], which is in 
general well-explained by the mechanism of curvature 
radiation. However, an extension of the spectral mea- 
surements for the brightest gamma-ray pulsars towards 
both, higher and lower energies, has revealed that the 
spectra beyond the cut-off could be smoother (b~0.5). 
For example, the phase- averaged spectrum of the Crab 
pulsar is better fitted with the combination of param- 
eters b = 0.43, T = 1.59 and E = 0.50 GeV [54], 
rather than b = 1, T = 1.97 and E = 5.8 GeV as 
reported earlier by the Fermi-LAT collaboration based 
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on smaller gamma-ray statistics [55]. In any case, if the 
above noted fit of the energy spectra is extrapolated to 
higher energies, a dramatic decrease of gamma-ray fluxes 
well beyond 10 GeV is expected, preventing the detec- 
tion of pulsed emission with the current instrument at 
^100 GeV. The MAGIC telescope, using a novel trigger 
system detected sub-lOOGeV pulsed emission from the 
Crab pulsar [56], favouring models with exponential or 
sub-exponential cut-offs (slot gap and outer gap models). 
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Fig. 4 Spectral energy distribution (SED) of the pulsed gamma- 
ray emission from the direction of the Crab pulsar and nebula. 
Fermi-LAT points are shown (blue squares) together with MAGIC 
(grey and pink) and VERITAS (green) points. A Fermi-L AT- 
points best-fit, using two different hypotheses (b = 1, Eq = 5.8 
GeV and T = 1.97 and b = 0.85, E = 7 GeV and T = 1.97), is 
displayed in grey. The pulsed VHE radiation can be successfully 
accounted for (light blue, blue, green and red curves) by inverse 
Compton up-scattering of the pulsed magnetospheric X-ray emis- 
sion by a cold ultra-relativistic pulsar wind (see Sec. 3.2). From 
Ref. [60]. 

Yet unexpectedly, pulsed 7-ray emission above 100 GeV 
and up to 400 GeV of unknown origin was recently de- 
tected from the Crab with the VERITAS and MAGIC 
telescopes [57,58], cf. Fig. 4, challenging models for the 
origin of the periodic emission in neutron stars. Differ- 
ent explanations could be pursued to accommodate these 
new experimental findings within current models, such 
as secondary emission of electrons in the outer magneto- 
sphere [59] or IC emission from energetic electrons in the 
ultra-relativistic pulsar wind [60] (cf. also Sec. 3.2). Ap- 
proaches like these predict different spectral shapes and 
light-curve behaviour at GeV and TeV energies. The de- 
tected phase- averaged, pulsed emission (Fig. 4) could in 
principle be fitted by extrapolating the reported Fermi 
fluxes to the VHE domain as a power law with photon 
index of 3.8 ± 0.5 and a flux of 1% of the flux of the 
Nebula at 150 GeV, but the nature of such an extrapo- 



lation seems rather difficult to justify on physical (mag- 
netospheric) grounds [60]. The VHE light curve shows a 
double peak structure well-aligned with the light curve 
at lower energies, although narrower by a factor of two or 
three than those measured by Fermi-LAT. The spectrum 
of the narrow peaks, extending no more that 10% of the 
rotational period, does not show a significant deviation 
in its shape from the global spectral fit. Assuming a 
common (magnetospheric?) origin, a smooth connection 
of the VHE points with the HE points can be achieved 
by fitting the data with a broken power-law function, but 
to the exclusion of an exponential cut-off. An alternative 
explanation consists in considering the entire gamma-ray 
region as a superposition of two separate components, a 
nominal (magnetospheric) GeV one and an additional 
VHE component produced by IC up-scattering of the 
magnetospheric emission by the fast pulsar wind [60]. 
Measuring the spectral shape with high precision in the 
near future will provide constrains on these models and 
allow to investigate the connection with the low-energy 
points around 50 GeV and the spectral extension above 
400 GeV. Up to now, the observed 7-ray features make 
the Crab a unique source of this kind at VHE. An in- 
crease of the sample by observing the brightest Fermi- 
LAT pulsars, such Vela or Geminga will be pursued by 
H.E.S.S. II, MAGIC II and VERITAS (and CTA in the 
future), providing more input to understand the origin 
of this pulsed VHE radiation [61]. 

2.3 Pulsar Wind Nebulae 

Relativistic winds from energetic pulsars carry most of 
the rotational power into the surrounding medium, ac- 
celerating particles to high energies, either during their 
expansion or at the shocks produced in collisions of the 
winds with the sub- sonic environment. Accelerated lep- 
tons can interact with magnetic fields and low-energy 
radiations fields of synchrotron, thermal or microwave- 
background origins. As a result, non-thermal radiation 
is produced from the lowest possible energies up to ^100 
TeV. For magnetic fields of few /iG, freshly injected elec- 
trons (and positrons) create a synchrotron nebula around 
the pulsar, ranging from the radio to the X-ray and, in 
some cases, to the MeV band. At high energies a sec- 
ond component appears as a result of Comptonization of 
these soft photon fields by the relativistic leptons, creat- 
ing an extended IC-nebula around the pulsar [56]. 
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Fig. 5 Spectral energy distribution (SED) of the Crab Nebula 
in the high- and very high energy gamma-ray domain. The spec- 
tral points from low to VHE gamma-rays are shown together with 
a fit of the synchrotron component (blue dashed line) and predic- 
tions for IC gamma-rays calculated for three different values of the 
mean magnetic field: B = 100 fiG (solid red line), B = 200 /iG, 
and the equipartition field of the nebula of 300 fiG. From Ref. [55], 
reproduced by permission of the A AS. 

VHE observations of these pulsar wind nebulae (PWNe) 
have revealed PWNe to be the most effective Galactic 
objects for the production of VHE gamma-rays, allow- 
ing the detection of such systems even outside our own 
Galaxy (in the LMC [63]). As recently as of 2004, only 
the Crab PWN was detected with a steady gamma-ray 
fluxabove ITeVof (2.1±0.1 st at)xl0- 11 cm- 2 s- 1 [64,65]. 
The development of the new sensitive IACTs in the last 
years has raised the number of likely PWNe detected 
to at least 27 sources, whereas many of the unidentified 
gamma-ray sources are widely believed to be PWNe (or 
old relic PWNe) [23]. 

For many years, the Crab nebula was considered as a 
standard candle for the cross-calibration of VHE detec- 
tors, as the brightest persistent point-like TeV gamma- 
ray source seen effectively from both hemispheres. The 
main features of its non-thermal emission, extending over 
21 decades of frequencies, has been satisfactorily de- 
scribed by the formation of a PWNe based, to a large 
extent, on a simple MHD model for the interaction of a 
cold ultra-relativistic electron-positron wind with the in- 
terstellar medium [66]. Recent detailed two-dimensional 
MHD simulations [67,68] have confirmed such a concept, 
at least for the Crab Nebula. The IC emission detected 
at TeV provides crucial information about the conditions 
in the nebula even when it only constitutes a small frac- 
tion of the synchrotron luminosity of the nebula. In par- 
ticular, a comparison of the X-ray and TeV gamma-ray 
fluxes observed from the Crab Nebula has lead to a ro- 
bust estimate of the average nebular magnetic field of 
less than 100 /iG, in good agreement with predictions for 
the termination of the wind in MHD theory [66]. Fig- 
ure 5 shows the high-energy coverage of the Crab Neb- 



ula spectrum. While the COMPTEL and EGRET data 
carry information about the synchrotron radiation in 
the cut-off region, the Fermi-LAT data reveal the sharp 
transition from the synchrotron to the IC component at 
around 1 GeV. At an energy E^lOOGeV, a clear indi- 
cation of the IC maximum is supported by both satel- 
lite (Fermi-LAT [55], and ground-based (MAGIC [69] 
and VERITAS [58]) measurements, which show remark- 
able agreement with each other. The measurements with 
ground-based IACTs have almost approached 100 TeV 
[64,70,71], where the IC component should still extend to 
the energy region set by the maximum energy of the ac- 
celerated electrons, i.e., 1 PeV. Although the production 
of gamma-rays at such energies takes place in the Klein- 
Nishina regime, and is therefore strongly suppressed, fu- 
ture instrument such CTA should be able to detect this 
emission. 

Yet, despite the large coverage and deep observations, 
many aspects of this unique source are still unresolved. 
For instance, rapid high-energy flares with rise time as 
short as 6 hours from the Crab PWN have been reported 
by the Fermi-LAT and the AGILE collaboration [54,72]. 
This amazing discovery has opened new questions such 
as how these flares connect with the pulsar energy re- 
lease or as to their origin (are they related to the in- 
ner pulsar wind or to the magnetosphere?, see e.g., [73- 
77]). The exceptionally high fluxes during the active 
state in April 2011 allow detailed spectroscopy for dif- 
ferent flux levels [54]. In order to study the spectral 
evolution of the flaring component, a steady- state (con- 
stant) background has been assumed with a steep power- 
law spectrum described by a photon index Tb = 3.9. 
The spectrum of the flaring component has been as- 
sumed in the form of power-law with exponential cutoff, 
vF v = foE 2 ~ ri exp[-(E/E ) K ]. The results show that 
the spectra during all selected windows can be well de- 
scribed by the same photon index Tf = 1.27 ± 0.12 and 
exponential cutoff index k = 1, but with variable total 
flux fo and the cut-off energy Eq. A variation by a fac- 
tor of two allows a good fitting of the data, but the total 
flux has to be changed more than an order of magni- 
tude in this approach. While different theories (includ- 
ing synchrotron radiation and reconnection) have been 
put forward to explain these flares, many key issues are 
still unresolved. 

Even as one of the strongest sources in the TeV sky, 
the Crab nebula is very inefficient in producing gamma- 
rays through IC scattering, and only its extremely high 
spin-down power compensates for this. The energy den- 
sity of the magnetic field (of the order of ~ 100 /iG) 
exceeds by more than two orders of magnitude the ra- 
diation energy density. Thus, less than one per cent of 
the energy of the accelerated electrons is released in IC 
gamma-rays, the rest being emitted through synchrotron 
radiation. In other systems, the pulsar wind is not as 
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powerful as the one in Crab, resulting in weaker mag- 
netic fields in the nebula of the order of a few fiG. This 
low magnetic field translates into a more efficient emis- 
sion via IC at VHE due to the sharing of the electron 
energy losses between synchrotron and IC mechanism. 
For instance, in the case of the cosmic microwave radi- 
ation (CMB), the two radiation components are related 
through L 7 /Lx = ^cmb/% — 1 (5/3/iG) -2 . This im- 
plies that in a PWN with a nebular magnetic field of 
about 10 fiG or less, the IC gamma-ray production ef- 
ficiency could be as large as 10%. Given that the ro- 
tational energy of pulsars is eventually released in rel- 
ativists electrons accelerated at the termination shock, 
PWNe associated with young pulsars with spin-down lu- 
minosities Lq ^ 10 34 (<i/lkpc) 2 erg/s were expected to 
be detected [78]. These expectations have been con- 
firmed by the results obtained with MAGIC and VERI- 
TAS, but overall by the survey performed with H.E.S.S. 
The Galactic plane survey (GPS) as seen by H.E.S.S. in 
Fall 2012 is shown in Fig. 6. The survey, covering a 
range between [-85°, 60°] in longitude and [-2.5°, 2.5°] 
in latitude, has revealed more than fifty new VHE 7-ray 
sources, out of which more than half are believed to be 
gamma-ray PWNe, located in the close vicinity of young 
and energetic pulsars. 



tracing the cooling mechanisms in the leptonic popula- 
tion injected by the pulsar (as observed in HESS J1825- 
137 or HESS J1303-631 [82,83], cf. Fig. 7). 
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Fig. 6 Significance (pre-trial) map of the Galactic plane survey 
by H.E.S.S. From Ref. [296]. 

Presently PWNe constitute the largest galactic TeV 
source population. Many previously dubbed "dark" TeV 
gamma-ray sources, including the first unidentified TeV 
gamma-ray source discovered by the HEGRA collabo- 
ration, TeVJ2032+4130 [79], have later been identified 
as PWNe. Most of these identifications with PWNe 
are quite convincing, yet still tentative, except for sev- 
eral ones which are firmly identified, either by excel- 
lent radio/X-ray morphological correlations, such as the 
Kookaburra complex, MSH 15-52 and VelaX [80,81], 
or by observations of an energy-dependent morphology, 




Fig. 7 VHE image of the TeV pulsar wind nebula candidate 
HESS 1303-631 at different energy ranges. The highest-energy 
photons originate near to the pulsar. X-ray (XMM) contours are 
shown in white. See Ref. [297] 

Out of the PWNe detected at VHE two different pop- 
ulations of PWNe seem to be emerging: PWNe associ- 
ated to young, compact X-ray PWNe, often still embed- 
ded in their associated supernova remnant; and evolved 
(extended and resolved) sources, in which the TeV emis- 
sion seems to be due to a "relic" population of electrons, 
whereas the associated shell has already faded away. In 
the latter group, the centre of gravity of the extended 
TeV images is often offset with respect to the position 
of the powering pulsar. Asymmetric, one-sided images 
of these PWNe have also been found in X-rays, but on 
significantly smaller scales. Although the mechanism 
which causes PWN offsets from the pulsar positions is 
not yet firmly established, this effect could be linked to 
the propagation of a reverse shock created at the ter- 
mination of the pulsar wind in a highly inhomogeneous 
medium [62]. The significantly larger extension of the 
TeV emission region can be understood as a result of 
several factors: (i) Generally, for PWNe with magnetic 
field of order of 10 /iG or less, as apparently the case for 
most TeV PWNe, the electrons responsible for the X- 
ray emission are more energetic than the electrons emit- 
ting TeV gamma-rays. Therefore, synchrotron-burning 
of the highest-energy electrons results in a smaller size 
of the X-ray source, (ii) When electrons diffuse beyond 
the PWN boundary, they emit less synchrotron radia- 
tion (due to the reduced magnetic field), but they can 
still effectively radiate gamma-rays via inverse Compton 
scattering of the universal CMB. (iii) Finally, because of 
the high X-ray background, the sensitivities of X-ray de- 
tectors like Chandra and XMM- Newton are dramatically 
reduced beyond several angular minutes. This signifi- 
cantly limits the potential of these instruments for weak, 
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extended X-ray sources. In contrast, the sensitivity of 
IACT arrays remains almost unchanged approximately 
within a 1° radius of field-of-view. This flat response 
makes IACT technique the most powerful tool for study- 
ing the non-thermal population of electrons in PWNe. 
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Fig. 8 Energy-dependent VHE morphology of pulsar wind neb- 
ula HESS J 1825-137, showing a softening of the spectra with in- 
creasing distance from the pulsar. The plot shows the energy spec- 
tra in radial bins as indicated in the inset (with the dashed line 
from the innermost region for comparison). From Ref. [82] 

The asymmetry observed in those PWNe has been ex- 
plained as a consequence of the propagation of the pre- 
cursor supernova explosion in the inhomogeneous inter- 
stellar medium [84] , resulting in a faster evolution of the 
associated PWN in the opposite direction of the denser 
environment or/and a high kick-off velocity of the pulsar, 
displacing it from the centre of the supernova explosion. 
The accumulation of particles with time, the continuous 
injection and the ubiquitous presence of a soft photon 
target (CMB) make these objects extremely efficient in 
the production of VHE emission. The high flux and ex- 
tension of these TeV PWNe have permitted the inves- 
tigation of the spectral behaviour with good statistics 
in different regions of the nebula, unveiling a softening 
of the gamma-ray spectral index as a function of the 
distance from the pulsar (see Fig. 8). This effect is 
due to the radiation of uncooled electrons which quickly 



leave the compact region near the pulsar, suffering signif- 
icant radiative losses as they propagate away. This seems 
also to be the case for Vela X, a nearby PWN related 
to the powerful pulsar PSR J0835-4510 (r ^11,000 yr, 
L = 7 x 10 36 erg/s). Vela X has been established [81] 
as one of the strongest TeV gamma-ray sources in the 
Galaxy. The energy spectrum of this source is quite 
different from other galactic sources; it is very hard at 
low energies, with photon index T & 1.5, and contains 
a high-energy exponential cut-off resulting in a distinct 
maximum in the SED at 10 TeV. Because of the nearby 
location of the source (d « 300 pc) we see, despite the 
large angular size of the gamma-ray image of order of 
1 degree, only the central region with a linear size less 
than several pc. In this regard, Vela X is a perfect object 
for the exploration of processes in the inner parts of the 
nebula close to the termination shock. The significantly 
improved sensitivity of the future CTA instrument and 
its superior angular resolution (one to two arc minutes at 
10 TeV) should allow a unique probe of the relativistic 
electrons inside the region of the termination shock, i.e., 
at the very heart of the accelerator. 
Along with these evolved nebula, a large number of 
compact objects have also been identified recently (see, 
e.g. [85,86]), in which the PWN is still expanding 
within the shell. A text-book example is the compos- 
ite SNRG327.1-1.1 (HESS J1554-550) [87], in which the 
detected TeV emission is spatially coincident with the X- 
ray and radio PWN, well inside the remnant. A similar 
case is the newly detected source HESS J1818-154 [88], 
embedded in the SNR G15.4+0.1. The latter was discov- 
ered after a long exposure of 145 h with a flux of 1.5% 
of the Crab Nebula flux, and no X-ray or radio PWN 
has been detected yet, allowing SNR G15. 4+0.1 to be 
identified as a composite SNR by means of VHE obser- 
vations only. Those objects display a very low magnetic 
field in comparison to the Crab Nebula of the order of 
a few /iG, compensating so the lower spin-down power 
luminosity with a particle-dominated wind, which allows 
an enhancement of the inverse- Compton emission at very 
high energy. 

2.4 TeV Binary Systems 

The number of TeV binary systems - sources emitting 
variable, modulated VHE emission composed of a mas- 
sive star and a compact object - has increased steadily 
in the last years, thanks to the large time coverage and 
the deep and uniform exposure of the Galactic plane by 
MAGIC, VERITAS and H.E.S.S. The TeV emission is 
believed to arise from the interactions between the two 
objects, either in an accretion-powered jet (microquasar 
scenario), or in the shock between a pulsar wind and a 
stellar wind (wind-wind scenario) (see e.g. [89-93], cf. 
also Sec. 3.2). In the microquasar scenario, particle 
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acceleration takes place in a jet which originates from 
an accretion disk. This scaled-down version of an ac- 
tive galactic nucleus opens the possibility to obtain sig- 
nificant insights into the mechanism of jet production. 
In the wind-wind scenario, on the other hand, parti- 
cle acceleration occurs in the interaction region between 
a ultra-relativistic pulsar wind and the dense radiation 
field provided by the companion star. Likewise, X-rays 
and high-energy components are expected due to radia- 
tive (synchrotron and inverse- Compt on) cooling of rel- 
ativists electrons accelerated at the termination shock 
[94,95]. 

Four periodic binary systems have been firmly identified 
at VHE (PSRB1259-63 [96], HESS J0632+057 [97,98], 
LS5039 [99] and LSI +61 303 [100-102]), whereas two 
more sources (HESS J1018-589 [103] and CygX-1 [104]) 
are less certain and still pending confirmation. The 
observed variability implies a compact emission region 
which translates into a point-like source morphology at 
a distance of 1 to 5 kpc. Indeed, the majority of point- 
like sources detected in the H.E.S.S. Galactic Survey 
have been identify as TeV binary systems. This uni- 
vocal identification is based on the observed VHE vari- 
ability/periodicity and correlations with flux variation 
at other wavelengths. They exhibit a maximum flux of 
~ 5 — 15% of the Crab Nebula flux and apparent simi- 
lar spectral indices (2.0 to 2.7), but the enlargement of 
the TeV (and GeV) binary sample has indicated a very 
diverse behaviour from one system to the other, demand- 
ing a detailed source-to-source investigation. 
The first TeV binary established was the pulsar-B2Ve 
star system associated to PSRB1259-63 (or LS 2833) in 
2004, which was anticipated before its detection [95]. In 
this system, a 48 ms pulsar is moving around a massive 
Be star, crossing its disk every 3.4 years, on a highly 
eccentric (e=0.87) orbit. The observations show a com- 
plex light curve, and the VHE emission can be satisfacto- 
rily explained in a pulsar-wind stellar-wind scenario, al- 
though the different year-to-year observations still chal- 
lenge current models. Moreover, the source exhibited a 
large post-periastron orphan flare at GeV energy that 
was not observed in the TeV range [105,106,290], which 
lasted approximately two weeks with an enhanced flux 
above 100 MeV at the level of 3 x 10" 10 erg cm~ 2 s _1 . 
Several scenarios have been proposed to account for 
this phenomenon, involving energy-dependent absorp- 
tion processes and/or Comptonization of the photon field 
provided by the star by the cold ultra-relativistic pulsar 
wind [107]. 

The second, very-long-period (~320 days)-system was 
discovery serendipitously in the H.E.S.S. survey, being 
one of the very few point-like (<2 ; ) unidentified sources. 
HESS J0632+057 was finally identified in a joint H.E.S.S. 
and VERITAS campaign, followed up by MAGIC ob- 
servations, through long-term X-ray observations with 



SWIFT, that succeeded in confirming its nature. Strong 
evidence for periodic X-ray variability with a very long 
period of 321+5 days has been reported [108], implying 
the discovery of a binary for the first time on the basis 
of TeV observations (see also Fig. 9). 
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Fig. 9 VHE observations of the binary system HESS J0632+057 
folded with a period of 321 days. The H.E.S.S. (circular markers) 
and VERITAS (open squares) measurements are shown in different 
colours for different observational periods. From Ref. [97]. 

The last two mentioned VHE binary systems, LSI 
+61 303 and LS 5039, show short-periodic orbital vari- 
ability, of the order of days, allowing a larger integra- 
tion of VHE data and deeper investigation of their light 
curve. However, they behave quite differently from each 
other. While LS 5039 (P~3.9 days) exhibits are very 
regular light curve, LSI+61 303, with a period of ~26.5 
days, shows a quite erratic behaviour, likely related with 
a 1667 super-orbital variability [109]. The nature of the 
compact object for both system is unknown: It could 
be anything from a 1.4 M Q neutron star to a (3.7)4 M Q 
black hole. No pulsation has been found in radio or X- 
ray searches. It seems likely, however, that any pulsed 
radiation would be absorbed in the optical-thick dense 
ambient due to Compton scattering [22]. 

These two binary systems have also been detected with 
the Fermi-LAT telescope above 100 MeV. The spectrum 
of LS 5039 shows a clear hardening in the 0.3 to 20 TeV 
region (see Fig. 10), while the GeV component shows a 
softening in inferior conjunction. On the other hand, at 
superior conjunction an opposite behaviour is observed. 
LSI +61 303 on the contrary, does not show variation of 
the spectral index, but its emission vanished after Oc- 
tober 2008, reappearing again in 2010, accompanied by 
a change in the high-energy flux with decrease of the 
orbital modulation in 2009 [111-113]. From the multi- 
wavelength data it is clear that more sophisticated sce- 
narios are needed to understand the acceleration and 
emission processes involved in these two sources. 

Finally, two more VHE regions have been associated 
with binary systems: MAGIC has reported a 4cr evidence 
for VHE emission from the direction of the Microquasar 
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Fig. 10 High-energy (Fermi-LAT) and VHE (MAGIC, VERITAS and H.E.S.S.) observations of LSI +61 303 (left) and LS 5039 (right), 
cf. Refs. [298,299]. The right figure shows the spectral data at inferior conjunction in red circles whereas the observations in the superior 
conjunction are shown in blue triangles. 



Cyg X-l [104], correlated with an increase in soft and 
hard X-rays, but this was not confirmed during later, 
similarly high X-ray flux flaring events; and the GeV 
16.5 days binary system 1FGL J1018. 6-5856 [114], coin- 
cident with the H.E.S.S. source HESS J1018-589. For 
the latter, no VHE variability has been discovered yet, 
making the association somewhat unclear. Deep obser- 
vations and uniform exposure in time with H.E.S.S. will 
help to clarify the origin of its VHE emission. 

2.5 Galactic Centre 

The Galactic Centre (GC) harbours many remarkable 
objects, including a few potential sites for particle ac- 
celeration and gamma-ray production, in particular the 
compact radio source Sgr A*, a suspected super-massive 
black hole located at the dynamical centre of the Galaxy. 
The GC contains a strong gamma-ray source (cf. Figs. 
11 and 12) with a broad-band spectrum that spans from 
100 MeV [115] to 30 TeV [116]. Assuming that gamma- 
rays from the entire interval are linked to the same 
source, the spectrum has an interesting form with several 
distinct features: Hard at low energies, with a photon in- 
dex T w 2.2, it becomes significantly steeper by Ar ^0.5 
above 2 GeV [115], then hardens again at TeV energies 
with a photon index r c± 2.1 and an apparent break or 
cut-off above 10 TeV (see Fig. 12). 




Fig. 11 The image of the several-hundred parsec region of the 
Galactic Centre in TeV gamma-rays (top: 7-ray count map; bot- 
tom: same map after subtraction of the two point sources). It 
contains a point like source (angular radius less than a few arc- 
minutes), the gravity centre of which coincides with an accuracy 
of 13 arc-seconds with the compact radio source Sgr A* (marked 
with black star) - a supermassive black hole at the dynamical centre 
of the Milky Way [120,121]. The second point-like source located 
about one degree away positionally coincides with the composite 
supernova remnants G09+0.1 [85]. A prominent feature of this re- 
gion is the ridge of diffuse emission tracing several well- identified 
giant molecular clouds (lower panel; cf. Ref. [122] for more details). 
This complex region contains some other, not yet firmly identified, 
"hot spots". 

Although the gamma-ray source spatially coincides with 
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the position of Sgr A* (see Fig. 11), the upper limit 
on the angular size of the TeV source of a few arc min- 
utes is still too large to exclude the link to other po- 
tential sources located within the central ^ 10 pc region. 
The detection of variability of the gamma-ray flux would 
greatly contribute to the localisation of the gamma-ray 
production region in Sgr A*. However, unlike the ob- 
servations at radio and X-ray wavelengths, no variability 
has been observed both at GeV and TeV energies. This 
disfavours, but still cannot discard Sgr A* as a possible 
gamma-ray source, especially given that several radiation 
mechanism, associated with the accretion flow, are capa- 
ble of explaining the reported gamma-ray fluxes [117]. 
Perhaps a more plausible site of gamma-ray production 
could be the central, dense extended region of radius 
of 10 pc. However, even in this scenario Sgr A* re- 
mains a potential source indirectly responsible for the 
gamma-ray signal through interactions of runaway par- 
ticles accelerated in Sgr A*, but later injected into the 
surrounding dense gas environment [118,119]. The anal- 
ysis of the combined Fermi-LAT and H.E.S.S. data show 
that the complex shape of the GeV- TeV radiation can 
be indeed naturally explained by the propagation ef- 
fects of protons interacting with the dense gas within 
the central 10 pc region [115,119]. A good agreement 
between the data and calculations is shown in Fig. 12, 
where the radial profile of the gas density has been care- 
fully taken into account. The flat spectra in the seg- 
ments of the proton spectrum around 1 GeV, and at 
TeV energies (below 10 TeV) have different explana- 
tions. While at GeV energies the protons are diffusively 
trapped, so that they lose a large fraction of their energy 
before they leave the dense 3 pc region, at TeV ener- 
gies they propagate rectilinearly. At intermediate ener- 
gies the protons start to effectively leave the inner 3 pc- 
region, and the steepening of the energy spectrum can be 
naturally referred to the energy- dependent diffusion co- 
efficient. What concerns the proton injection spectrum, 
it should be a hard power-law, close to E~ 2 , with an in- 
trinsic cut-off around 100 TeV. The required total energy 
of protons currently trapped in the gamma-ray produc- 



tion region, W p 



lO^n/lO^cm 3 )" 1 erg 



is quite modest, given that the density in the circum- 
nuclear ring could be as large as 10 5 cm -3 [119]. 
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Fig. 12 Energy spectra of gamma-ray emission from GC. The 
Fermi-LAT [115] and H.E.S.S. data [116] are shown together with 
calculations of 7-rays from pp-interactions within radial cones of 
various size up to 50 pc. The flux falls off rapidly after 3 pc because 
the main contribution comes from the 1.2-3 pc circum-nuclear ring. 
From Ref. [119], reproduced by permission of the AAS. 

The interpretation of the spatially unresolved gamma- 
ray emission towards Sgr A* by interactions of runaway 
protons with the dense gas in the central (several pc) 
ring, predicts a smooth transition to another radiation 
component formed in more extended regions of the GC. 
The energy and spatial distributions of this radiation de- 
pend on the injection history of protons and the charac- 
ter of their diffusion. The H.E.S.S. observations of the so- 
called Central Molecular Zone (CMZ) of radius « 200pc 
indeed revealed an extended TeV gamma-ray emission 
[122] with a clear correlation with the most prominent 
giant molecular clouds located in CMZ (see Fig. 11). 
Using the maps of TeV gamma-ray emission, and maps 
of the CS (J=l-0) emission which contain information 
about the column density in dense cores of molecular 
clouds, the cosmic-ray density in these clouds has been 
derived. It appears to be significantly enhanced (by an 
order of magnitude at multi-TeV energies) relative to the 
local cosmic-ray flux in the solar neighbourhood. This 
indicates to a strong non-thermal activity accompanied 
with proton acceleration which in the past was perhaps 
higher than at the present epoch. An additional support 
for this hypothesis comes from the spatial distribution 
of gamma-rays. The H.E.S.S. observations show that 
the ratio of gamma-ray flux to the molecular gas column 
density varies with galactic longitude, with a noticeable 
"deficit" at I « 1.3°. This interesting feature can be in- 
terpreted as a non-uniform spatial distribution of cosmic 
rays, i.e. the relativistic protons accelerated in Sgr A* 
have not yet had time to diffuse out to the periphery of 
the 200 pc region. The epoch of the high activity of the 
accelerator depends on the proton diffusion coefficient. 
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Assuming, for example, that the propagation of multi- 
TeV protons in the GC proceeds with a speed similar to 
the one in the Galactic Disk, the epoch of high activity of 
the accelerator and the total energy release in relativistic 
particles during the outburst are estimated to be 10 4 yr 
and 10 50 erg, respectively [122]. 

High-energy processes that take place in the GC ap- 
parently play a key role in the formation of two enor- 
mous gamma-ray structures recently discovered in the 
Fermi-LAT data set - the Fermi bubbles [123]. Centred 
on the core of the Galaxy, these structures symmetri- 
cally extend to approximately 10 kpc above the Galac- 
tic plane. The parent relativistic particles (e.g., pro- 
tons) could be accelerated in the nucleus of GC, and 
then injected into Fermi Bubbles. Alternatively, protons 
and electrons could be produced in situ through first- 
and/or second-order Fermi acceleration mechanisms sup- 
ported by hydrodynamical shocks or plasma waves in 
a highly turbulent medium. The processes that create 
and support these structures could originate either from 
an AGN-type activity related to the central black hole 
(Sgr A*) or from ongoing star formation in the galactic 
nucleus. 

The luminosity of gamma-rays with hard, £? -2 -type, 
spectrum in the energy interval 1-100 GeV (see Fig. 13) 
is L 7 « 4 x 10 3T erg/s. Given the overall limited energy 
budget of the GC, particle acceleration and gamma-ray 
emission in the Fermi bubbles should proceed with very 
high efficiency. Despite the significant differences of the 
models proposed for the origin of the Fermi bubbles, only 
two radiation mechanism can be responsible for gamma- 
rays - IC emission by relativistic electrons or decays of 
neutral pions produced in pp-inter actions. Because of se- 
vere radiative energy losses, however, the mean free path 
of ^ 100 GeV electrons is significantly shorter than the 
size of the Fermi bubbles. Therefore, one has to pos- 
tulate in situ electron acceleration throughout the entire 
volume of the bubbles [123,124]. Such a scenario could be 
realised through stochastic (second-order Fermi) acceler- 
ation [125] or due to series of shocks propagating through 
the bubbles and accelerating relativistic electrons [126]. 
Importantly, the suggested acceleration mechanism seem 
unable to boost the electron energy beyond 1 TeV, thus 
in order to explain the extension of the observed gamma- 
ray spectrum up to 100 GeV by IC, one has to invoke 
FIR and optical/UV background emission supplied by 
the galactic disk (see Fig. 13). This model provides ro- 
bust predictions. In particular, since the FIR and opti- 
cal/UV contributions to the target field for IC scattering 
decrease quickly with distance from the disk, the spec- 
trum of gamma-rays from high latitudes should contain a 
cut-off above tens of GeV. The limb brightening at high- 
est energies is another characteristic feature predicted by 
this model. These spectral and spatial features can be 
explored in the near future, after the gamma-ray photon 



statistics in the Fermi-LAT data set has achieved an ad- 
equate level. 

An hadronic origin for the observed gamma-rays is an 
alternative interpretation suggested for the Fermi bub- 
bles [124,127]. Despite the low plasma density in the 
Fermi bubbles, n ^ 10 -2 cm -3 , the efficiency of pro- 
ton interactions can be very high. Indeed, if protons 
would have been continuously injected and trapped in 
the bubbles over timescales of approximately 10 10 yr, 
the main power in accelerated protons would be lost in 
pp- collisions given that the characteristic time of the lat- 
ter, t pp = l/(k p na pp c) « 5 x 10 9 (n/10" 2 cm" 3 ) -1 yr, 
is shorter than the confinement time. This implies that 
one deals with a so-called "thick target" scenario, when 
the system is in saturation. The hadronic gamma-ray 
luminosity is equal to L 1 « W p /t pp ^o 1 where W p is the 
total energy of protons in the bubbles, and t 



PP— ^7T l 



is the 

timescale for neutral pion production in pp-inter act ions. 
In the saturation regime, W p = Q p t pp (with Q p the in- 
jection rate of protons), assuming that the energy dissi- 
pation through pp-collisions is the dominant loss process. 
Since t 



pp 



1/3 tpp^^o, we have L 7 



Qp/3, thus about 
a third of the power injected into relativistic CRs emerges 
in gamma-rays (of all energies) independent of the local 
density, interaction volume and the injection time. Note 
that since the timescale of pp-interactions is comparable 
to the supposed age of the bubbles of 10 10 yr, the effi- 
ciency would be somewhat less. Also, one should take 
into account that at low energies, ionisation and adia- 
batic losses of protons play a non- negligible role, thus 
the overall efficiency for a broad energy spectrum of pro- 
tons would be reduced to several percent. The fluxes 
of hadronic gamma-rays shown in Fig. 13 confirm these 
simple estimates. Note that independent of the history of 
injection of relativistic protons, the current total energy 
in protons should be as high as W p = L^pp^o ~ 10 55 
erg which is comparable to the magnetic field energy in 
the bubbles (cf. Ref. [301]). 
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Fig. 13 The spectral energy distribution of gamma-rays from the 
Fermi bubbles compared to theoretical predictions, (i) IC model 
of Ref. [125] (solid line) assuming stochastic acceleration of elec- 
trons in the bubbles (the contributions from the scattering on the 
CMB, FIR, and optical/UV backgrounds are shown separately); 
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(ii) IC model of Ref. [126] (dotted line) assuming diffusive shock 
acceleration of electrons; (Hi) hadronic model of Ref . [124] (dashed 
line). The figure is from Ref. [125]. 

The above noted hadronic model of gamma-ray emis- 
sion of the Fermi bubbles does not exclude other 
"hadronic" scenarios with faster energy release related 
to the activity of the central black hole Sgr A*. A fast 
energy release can be provided, for example, by the cap- 
ture of stars by Sgr A* over the last 10 Myr with an 
average capture rate of3xl0 _5 yr _1 and energy release 
of 3 x 10 52 erg per capture [128]. It has been argued 
in Ref. [140] that quasi-periodic injection of hot plasma 
could produce a series of strong shocks in the Fermi bub- 
bles which can (re) accelerate protons beyond the "knee" , 
up to energies of about 10 18 eV. If confirmed by indepen- 
dent detailed hydrodynamical simulations, this could ap- 
pear a viable solution for the origin of one of the most 
"problematic" (poorly understood) energy intervals of 
cosmic rays. 

2.6 Blazars 

Most of the detected extragalactic gamma-ray sources 
belong to the blazar class, which comprises BL Lac ob- 
jects and Flat Spectrum Radio Quasars (FSRQs). The 
central engine in these active galaxies (AGNs), a super- 
massive black hole (BH) of mass ^ 1O T M0 surrounded 
by an accretion disk, is commonly believed to eject a 
relativist ic jet pointing almost directly towards the ob- 
server. Doppler boosting effects results in strong flux 
amplification, thus naturally favouring the detection of 
blazars on the extragalactic sky. 

Fermi-L AT, for example, has detected over 1000 extra- 
galactic high-energy (HE) sources in two years of sur- 
vey (2LAC), most (> 90%) of which are blazars [129]. 
In comparison, non-blazar sources like starburst galax- 
ies (SBs) or radio galaxies (RGs) only make out a minor 
fraction (in numbers). 

At the time of writing more than 50 extragalactic VHE 
sources, populating the whole sky, are listed in the 
online TeV Catalog (TeVCat).^ The majority of them 
(~ 90%) are again of the blazar type, with the so-called 
high-frequency-peaked BL Lac objects (HBLs, with low- 
energy component peaking at v v > 10 15 Hz, in con- 
trast to LBL s= low- frequency-peaked BL Lacs, peaking 
at v v < 10 14 Hz) constituting the dominant (^ 70%) 
sub-class, yet also including three FSRQs (3C279 at 
z = 0.536; PKS 1510-089 at z = 0.361 and PKS 1222+21 
at z = 0.432). FSRQs are typically distinguished from 
BL Lac objects by the presence of strong and broad (rest- 
frame equivalent width ^ 5 A) optical emission lines. Al- 
most all Fermi-detected FSRQs for which v v can be es- 
timated are of the low-frequency-peaked (v p < 10 14 Hz) 



type. Note that AGNs, which have been detected at TeV 
are typically characterised by a harder GeV photon in- 
dex than the majority of 2LAC sources. 
At present, blazar sources out to redshift z ~ 0.6 (i.e., 
3C279 at z = 0.536 [163] and BL Lac KUV 00311-1938 
at z > 0.51, tentative z = 0.61 [164]) have been detected 
at VHE energies, cf. Fig. 14 for their redshift distribu- 
tion. Blazar population studies at lower (radio-X-ray) 
frequencies indicate a redshift distribution for BL Lacs 
objects that seems to peak at z ~ 0.3, with only few 
sources beyond z ~ 0.8 (under the proviso of some bias 
as for a substantial fraction of BL Lacs the redshift is not 
known), while the FSRQ population is characterised by 
a rather broad maximum between z ~ (0.6 — 1.5) [160]. 
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Fig. 14 Distribution of redshift for the VHE-detected blazars. 
Redshift data are taken from TeVCat. Most objects are within 
z < 0.2. 

Observed VHE flux levels typically range from ~ 1% 
(for the average/steady states) up to ~ 15 times the 
Crab Nebula steady flux in high activity stages (e.g., for 
the 2006 flaring state of PKS 2155-304 [137]). Gamma- 
ray emission beyond 10 TeV has been established, with 
measured photon energies reaching, e.g. ~ 20 TeV in 
Mkn 501 (z=0.034; 1997 VHE flare [130]) and Mkn 421 
(z=0.031; 2010 VHE flare [131]). The observed VHE 
spectra usually vary between hard and soft power laws; 
for the HBL sources, for example, inferred photon indices 
(assuming a single power law dN/dE oc E~ r ) range from 
~ 2.3 to ~ 4.5, with some indications for spectral hard- 
ening with increasing activity. Non-HBL sources (i.e., 
IRL=intermediate-peaked BL Lacs and LBLs) are usu- 
ally detected during high states only, with low states 
expected to fall below current sensitivities. 
In a vF v representation (SED), blazars often show a 
double-humped structure, cf. Figs. 15 and 16. The first 
hump (sometimes reaching peak frequencies up to ~ 10 19 
Hz, as in the HBL object RGB J0710+591 [132]) is com- 
monly attributed to non-thermal electron synchrotron 
emission; the second one is widely believed to be due to 
leptonic inverse Compton processes (on the synchrotron 
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photons in SSC, or on ambient photons in External In- 
verse Compton [=EC] models), although hadronic sce- 
narios often remain possible. Different blazar population 
studies seem to suggest that there is a continuous spec- 
tral trend (see Fig. 15) from FSRQ -^LBL^IBL^HBL, 
often called the " blazar sequence" , characterised by a de- 
creasing source luminosity, increasing synchrotron peak 
frequency and a decreasing ratio of high- to low-energy 
component [133,134] (but cf. also [135] for caveats due 
to selection effects and unknown redshift). 



extremely broad (stretching from keV to TeV), defying 
a simple homogeneous SSC interpretation [145]. 
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Fig. 15 Sequence of characteristic blazar SEDs as a function of 
source luminosity from FSRQ (top curve) to HBL objects (bottom 
curve). From Ref. [134], reproduced with permission ©ESO. 

Blazar SEDs can span almost 20 orders of magni- 
tude in energy, making simultaneous mult i- wavelength 
observations a particular important diagnostic tool to 
disentangle the underlying non-thermal processes. A 
variety of leptonic and hadronic emission models have 
been discussed in the literature (see, e.g., [136] and refer- 
ence therein). A significant correlation between TeV and 
X-ray flux variations for example, which is often found, 
could favour a leptonic synchrotron- Compton interpre- 
tation, but counterexamples ("orphan TeV flares") do 
exist [141]. Short-term variability is usually more diffi- 
cult to account for in hadronic models because of longer 
cooling timescales, but strong magnetic fields (for proton 
synchrotron, e.g. [142]) or high target matter densities 
(pp-interactions triggered by jet-star interactions, e.g. 
[143]) may partly compensate. While for HBL objects, 
homogeneous (one-zone) leptonic SSC modeling often 
seems to provide a reasonable SED characterization (but 
see, e.g., [144] for a possible exemption), this does not 
apply in a similar way to LBL objects. Among the four 
LBLs detected, for example, AP Lib (z = 0049) repre- 
sents an intriguing example where the 2nd bump seems 
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Fig. 16 A recent, double- hump-structured SED example: The 
high-frequency-peaked (HBL) BL Lac object PG 1553+113 as 
based on VHE (MAGIC, 2005-2009) observations and archival 
data. Pronounced variability (on yearly time scale) is seen in the 
X-ray band. The average SED has been modelled with a one-zone 
SSC model (continuous black line). From Ref. [159], reproduced 
by permission of the A AS. 
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FSRQ 
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~ 10 min 
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Tab. 2 VHE variability in AGN: Characteristic minimum VHE 
variability timescale A^vhe as observed with current instruments 
for an exemplary set of AGN. 

Despite the limited temporal coverage of the current 
IACTs more than half of the AGN detected in the TeV 
domain shown variability, albeit often weak. For the ma- 
jority of them, variability timescales above one month 
have been found. In about a quarter of them there 
is clear evidence for short-term VHE variability on ob- 
served timescales of less than one day, cf. Table 2. The 
HBL class currently reveals the most rapid and dramatic 
VHE gamma-ray flux variability with observed variabil- 
ity timescales < 5 min, as found by the H.E.S.S. and 
MAGIC experiments for PKS 2155-304 (z = 0.116) [137] 
and Mkn 501 [138], respectively, cf. Fig. 17. Given the 
limited angular resolution (~ 0.1°) of IACTs, this im- 
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plies that one of the most constraining requirements on 
the jet kinematics and the high-energy emitting region 
comes from VHE variability studies. Fast VHE variabil- 
ity from distant blazars can also be used to derive con- 
straints on an energy-dependent violation of Lorentz in- 
variance (energy-dependent speed-of-light) as predicted 
in various models of Quantum Gravity [146,158]. 




Time - MJD53944.0 [min] 

Fig. 17 Light curve: Integrated flux /(> 200 GeV) versus time 
as observed by H.E.S.S. for PKS 2155-304 on July 28, 2006. The 
data are binned in 1-minute intervals. The horizontal line gives the 
steady flux from the Crab Nebula for comparison. From Ref. [137]. 

The detection of a large number of gamma-ray emit- 
ting blazars has opened a new research area - "obser- 
vational gamma-ray cosmology". The underlying idea 
is based on the energy-dependent absorption of 7-rays 
from distant extragalactic objects caused by interactions 
(7vhe 7ebl — > e + e~) with the Extragalactic Back- 
ground Light (EBL) that extends from UV to far IR 
wavelengths. The identification of absorption features in 
the spectra of 7-rays above 10 GeV, as well as detec- 
tion of characteristic angular and time distributions of 
gamma-rays produced during the cascade development 
in the intergalactic medium on large (^100 Mpc) scales, 
should allow us to derive unique cosmological informa- 
tion about the EBL and the intergalactic magnetic fields 
(IMFs). The realization of these exciting possibilities re- 
quires not only precise spectroscopic measurements from 
a large number of extragalactic objects located at dif- 
ferent redshifts, but, more importantly, a good under- 
standing of the intrinsic gamma-ray spectra. So, far 
the most significant contribution in this area comes from 
the measurements of gamma-rays from blazars with red- 
shifts between 0.1-0.2. In particular, based on such ob- 
servations, the H.E.S.S. collaboration has first reported 
a quite meaningful upper limit on the EBL at near 
and mid-infrared wavelengths [147]. Remarkably, the 
inferred upper limit appeared to be very close to the 
lower limit given by the measured integrated light of re- 
solved galaxies (galaxy counts), cf. also [148,150,151] 
for related inferences. Very recently, a similar result has 
been reported by the Fermi-L AT collaboration [152] for 
the EBL at optical and UV bands. One should men- 



tion, however, that the inferred upper limits are not 
model-independent. The H.E.S.S. result, for example, 
is based on the assumption that the differential intrin- 
sic spectrum is not harder than E~ 1,5 . The Fermi-L AT 
result is based on the detection of cutoffs in the aver- 
aged spectra of three samples of BL Lac objects com- 
bined in three different intervals of redshift, assuming 
that these cut-offs are caused by intergalactic absorption. 
Although both assumptions sound quite reasonable, and 
the derived upper limits agree with most of the theo- 
retical/phenomenological predictions for the EBL, one 
should keep in mind that they are not free of model as- 
sumptions. It is believed that future measurements by 
next-generation detectors, in particular by CTA, based 
on a much larger sample of AGN should significantly 
increase the source statistics and improve the quality 
of data, and consequently reveal details in the EBL. 
This optimistic view may, however, underestimate the 
difficulties related to the uncertainties of the intrinsic 
source spectra. On the other hand, the limits in which 
presently the EBL fluxes are robustly constrained, are 
quite tight, so one can "recover" the intrinsic gamma- 
ray spectra with a reasonable accuracy. Interestingly, in 
the case of some blazars, the gamma-ray spectra after 
correction for intergalactic absorption, appear extremely 
hard with photon indices ^ 1.5 or even close to 1, see 
e.g. [147,150,153]. This challenges conventional radia- 
tion models, but still cannot be considered as a failure 
of the standard blazar paradigm and as a need for new 
physics. Such spectra can still be explained, assuming, 
for example, the prevalence of certain conditions for the 
formation of the parent electron spectra (e.g., stochastic 
acceleration) or specific internal gamma-ray absorption, 
see e.g. [154-156]. Nevertheless, the growing number of 
VHE blazars with redshift exceeding z ~ 0.5 tells us that 
one should perhaps be prepared for even more dramatic 
assumptions, including violation of Lorentz invariance 
or "exotic" interactions involving hypothetical axion-like 
particles. An alternative interpretation of gamma-rays 
from very distant blazars (in case of their detection) ex- 
ists in the framework of standard physics: TeV gamma- 
rays can in principle be observed even from a source at 
z ^ 1, if the observed gamma-rays are secondary pho- 
tons produced in hadronic interactions (with CMB or 
EBL background photons) of energetic cosmic-ray pro- 
tons, originating in the blazar jet and propagating over 
cosmological distances almost rectilinearly. In the case 
of a detection of TeV gamma-rays from a blazar with 
z ^ 1, this model could in principle provide a viable 
interpretation consistent with conventional physics, but 
with an extreme assumption on the strength of the IMF 
in the range of 10~ 17 - 10~ 15 G (see, e.g. [157]). On the 
other hand, if VHE 7-rays from distant blazar attenuate 
through pair-production with EBL photons, constraints 
on the strength of the IMF can be derived by model- 
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ing the anticipated GeV emission from the electromag- 
netic cascades, taking the possible deflections of pairs in 
the IMF into account. According to a recent study, this 
method suggests a lower bound on the IMF of B ^ 10 _1T 
G [149]. 

2.7 Radio Galaxies 

Misaligned (non-blazar) AGNs, characterised by jets 
substantially inclined with respect to the observer, rep- 
resent a particularly interesting class of VHE emitters. 
Nearby radio galaxies (RGs) are especially attractive as 
their proximity may allow us to resolve the radio jets 
down to sub-parsec scales and to study possible multi- 
wavelength correlations. The absence of strong Doppler 
boosting could make a VHE detection challenging, yet 
also allow to get unique insights into emission regions 
otherwise hidden. 

Out of ~ 1000 high-energy (HE) sources (886 in the 
"Clean Sample"), Fermi-L AT has reported the detec- 
tion of only about ten misaligned RGs at GeV energies, 
with a predominance of the Fanaroff-Riley-type I (FR I) 
[129,161,162]. At TeV energies, only four RGs have been 
identified by current IACTs: The nearest AGN Cen A 
(d ~ 3.8 Mpc), the giant RG M87 (~ 16.7 Mpc), and the 
Perseus Cluster (d ~ 77 Mpc, z ~ 0.018) RGs NGC 1275 
and IC 310. A detection of the RG 3C66B was initially 
reported by MAGIC (2007 observations [139]), but the 
VHE emission seems not sufficiently disentangled from 
the nearby (separation ~ 0.12°) IBL blazar 3C66A to 
include it here. 

Cen A was detected at VHE in a deep (>120h) expo- 
sure by H.E.S.S. with a integral flux above 250 GeV of 
~ 0.8% of the steady flux of the Crab Nebula (corre- 
sponding to an apparent isotropic luminosity of L(> 250 
GeV) ~ 2 x 10 39 erg/s) [165]. The measured VHE spec- 
trum extends up to ~ 5 TeV and is consistent with a 
power-law of photon index 2.7 ±0.5. No significant vari- 
ability has been found. Fermi-LAT has also detected 
HE emission up to 10 GeV from the core of Cen A, with 
the HE light curve (15 d bins) being consistent with no 
variability and the HE spectrum described by a com- 
parable photon index [166]. A simple extrapolation of 
the Fermi HE power-law to the VHE domain, however, 
tends to under-predict the observed TeV flux. This could 
be indicative of an additional contribution to the VHE 
domain beyond the common synchrotron-Compton emis- 
sion, emerging at the highest energies [167,303]. While 
the giant radio lobes are also detected at GeV energies 
(with evidence for a spatial extension beyond the radio 
image [168]), they are clearly excluded (given the angu- 
lar resolution of H.E.S.S.) as source of the detected TeV 
emission. 

The giant radio galaxy M#7 was the first RG detected at 
TeV energies [169]. Commonly considered as a FR I- type 



RG, M87 is known to host a highly massive black hole of 
^bh — (2 — 6) x 10 9 Mq and to exhibit a relativistic jet 
misaligned by an angle ~ (15 — 25)°, consistent with 
modest Doppler boosting D = l/[Fj(l — j3 cos 9)] ^3 
(for a recent review, see e.g. [170]). M87 has shown par- 
ticularly interesting VHE features over the past years, 
including rapid variability (day-scale flux doubling time 
scales At bs ~ 1) above 350 GeV during high source 
states (with levels sometimes exceeding 10% of the Crab 
Nebula) , and a hard spectrum (compatible with a power 
law of photons index Y ~ 2.2) extending beyond 10 TeV 
[171-174]. Both the hard VHE spectrum and the ob- 
served rapid VHE variability are remarkable features for 
a misaligned AGN. M87 is the only RG where clear evi- 
dence for day-scale TeV variability (see Fig. 18) has been 
found, although there are now hints for day-scale vari- 
ability in IC 310 as well. The VHE variations in M87 are 
the fastest observed in any other waveband so far, and 
already imply a compact size of the 7-ray emitting region 
(R < Z)cAt bs/[l + 2])? comparable to the Schwarzschild 
radius (r s = (0.6 — 1.8) x 10 15 cm) of the black hole in 
M87. Consistent with radio evidence for the ejection of 
a new component from the core during a flare in 2008 
[175], this could point to the black hole vicinity as the 
most likely origin of the variable VHE radiation. 
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Fig. 18 VHE flare of the radio galaxy M87 as observed in April 
2010 by different IACTs. Significant day-scale activity is evident. 
The curve shows a fit of an exponential function to the data. From 
Ref. [174]. 

IC 310 and NGC 1275 are relatively new additions, be- 
ing discovered at VHE as part of a Perseus galaxy cluster 
campaign by MAGIC [176,178] . IC 310, has been origi- 
nally classified as a head-tail RG. However, according to 
a recent high-resolution VLBA study, there are little in- 
dications for jet bending [180], rather supporting (along 
with the jet-to-counter-jet flux ratio) the case that IC 310 
may instead rather be a weakly beamed blazar [177] . The 
source has been detected above 300 GeV in about 21h of 
data (taken in 2009/2010) at an average level of ~ 3% of 
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the Crab Nebula [176]. The VHE spectrum between 150 
GeV and 7 TeV is very hard (even harder than in M87) 
and compatible with a single power law of photon index 
T ~ 2.0. There is clear evidence for VHE variability on 
yearly and monthly time scales, with indications for day- 
scale activity found in a new analysis, features that are 
all reminiscent of the VHE activity seen in M87. 
On the other hand, the central dominant (FR I) clus- 
ter galaxy NGC 1275 (having radio jets misaligned by 
£ 30°), has been recently detected above ~ 100 GeV 
during enhanced high energy (Fermi-LAT) activity in 
46h of data (taken between 08/2010-02/2011). While the 
Fermi-LAT data reveal evidence for flaring activity above 
0.8 GeV down to time scales of days [179], the situation 
at VHE energies is less evident. No evidence of variabil- 
ity has been found in the 08/2010 to 02/2011 VHE light 
curve. A recent, improved analysis of an earlier (10/2009 
-02/2010) data set, however, seems to provide hints for a 
possible month-type VHE variability. NGC 1275 shows 
a steep VHE spectrum (r ~ 4.1) extending up to ~ 500 
GeV [178] and a hard HE (Fermi-LAT) spectrum (pho- 
ton index T ~ 2.1), indicative of a break or cut-off in the 
SED around some tens of GeV. 

2.8 Starburst Galaxies 

Starburst Galaxies (SGs) are galaxies showing a very 
high rate of star formation ("starburst") in a localised 
region, the burst sometimes being triggered by a close 
encounter with another galaxy. The resultant highly 
increased supernova (SN) explosion rate and the ex- 
pectation that the remnants (SNR) of those are effi- 
cient cosmic-ray (CR) proton accelerators (possibly up 
to ~ 10 16 eV [184]), suggest that starburst regions may 
possess a high cosmic-ray density. Because of the very 
high ambient gas densities (n > 100 cm -3 ), hadronic 
interactions (inelastic proton-proton collisions and sub- 
sequent 7r°-decay) could then lead to efficient 7-ray pro- 
duction, making SGs promising targets for HE and VHE 
astronomy. 

The spiral galaxy NGC 253 is the closest (d ~ 2.6 - 3.9 
Mpc) SG in the southern sky, harbouring a compact 
(~ 100 pc in extension) starburst region with an esti- 
mated SN rate of ~ 0.03 year -1 . The mean density of 
interstellar gas in this region (n ~ 600 protons/cm 3 ) is 
almost three order of magnitude higher than the average 
gas density in the Milky Way. NGC 253 was detected 
above 220 GeV in a deep exposure (> 119h) in 2005- 
2008 by H.E.S.S. at a flux level of - 0.3% of the Crab 
Nebulae [181]. A recent spectral analysis of an extended 
(-2009) data set (Fig. 19) shows that the VHE spec- 
trum can be described by a power law with photon index 
T c± 2.2, being compatible with a power-law extension of 
the HE (Fermi-LAT) 7-ray spectrum and little evidence 
for a spectral break or turnover between 200 MeV and 



5 TeV [182]. 
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Fig. 19 Differential energy spectrum of the starburst galaxy 
NGC 253 from high-energy (Fermi-LAT, crosses) to the VHE 
(H.E.S.S., circles). The shaded area shows the la confidence band 
from a simultaneous fit to the high-energy and VHE data. From 
Ref. [182]. 

M82, on the other hand, is a well-known SG in the North- 
ern Sky at a distance of d ~ 3.9 Mpc, harbouring an 
active starburst region with a high SN rate of 0.1 — 0.3 
year -1 . This source was detected above a 700-GeV anal- 
ysis threshold in a deep exposure (> 137 h) in 2008-2009 
by VERITAS at a level of 0.9% of the Crab Nebula (cor- 
responding to L 7 (> 700 GeV) ~ 2 x 10 39 erg/s) [183]. 
The noted SG findings are roughly compatible with 
hadronic CR (pp-interaction) models and energy- 
independent (advective) particle transport from the star- 
burst region (but see also [304,305] for an alternative 
approach related to the TeV emission from PWNe). If 
CRs would lose most of their energy due to inelastic 
collisions before they escape from the starburst region 
("quasi-calorimeter" scenario), the energy required in 
TeV protons would be of the order W p ~ L^tpp ~ 
10 52 (200cm _3 /n) erg, with L 1 ~ 10 39 erg/s the char- 
acteristic VHE 7-ray luminosity, t pp ~ (ccippkpn) -1 
the typical pp-collision time scale, and n the ambient 
gas number density. This would imply a typical CR 
energy density W p /V ~ 10 -10 erg/cm 3 , i.e., roughly 
two orders of magnitude higher than the average value 
(~ leV/cm 3 ) for the Milky Way. Note however, that in 
starburst region is likely to be mildly calorimetric, with 
only ~ 20 — 30% of the CR energy channeled into pion- 
production [182]. 



2.9. Candidates: 
Black Holes 



GRBs, Galaxy Clusters and Passive 



Several additional source classes are thought to be 
possible VHE emitters, but have not yet been detected. 
This includes GRBs, Galaxy Clusters and Passive Black 
Holes: 

Gamma- Ray Bursts (GRBs) have been observed to emit 
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intense flashes of 7-rays up to a least some tens of GeV 
energies, and their intrinsic photon spectra are likely 
to extend to TeV energies due to synchrotron-inverse 
Compton processes (by electrons accelerated at internal 
or external shocks) in leptonic models or electromag- 
netic cascades (triggered by py- and pp-interactions from 
CR-accelerated protons) in hadronic models [185]. Ex- 
tending the energy range for GRBs to the VHE domain 
could thus be particularly important to disentangle the 
origin of their high-energy emission. Significant 7-ray 
emission is expected to occur on time scales from second 
to at most hours, so that VHE instruments need to be 
re-posited very quickly to catch the dominant part of 
the outburst. The VHE prospects are, yet, generally 
constrained by (i) 77-pair production within the source 
(e.g., see the Fermi-L AT evidence for a spectral break 
around one GeV in GRB 090926 A [189]), and (ii) in 
particular, interactions with EBL photons, limiting the 
detectability of VHE GRBs from cosmological distances 
to relatively nearby sources [z ^ 0.6). Despite several at- 
tempts, no significant VHE excess from GRBs has been 
found yet and only upper limits have been reported (e.g., 
[186-188]). It will be highly interesting to see whether a 
VHE detection becomes possible in the near future with 
the new HESS II instrument, having a drive system de- 
signed as to reduce the response time to GRB alerts and 
a lower energy threshold as to reduce EBL absorption. 

Massive Cluster of Galaxies, harbouring powerful AGNs 
and being surrounded by strong accretion shocks repre- 
sent another potential VHE source class. Gamma-ray 
emission could be produced by several processes [190], 
e.g., IC up-scattering of CMB and other soft photon 
fields by electrons injected into the inter-cluster medium 
and further accelerated by turbulent processes [191,192]. 
Alternatively, hadronic (pp) interactions of relativistic 
CR protons with the ambient cluster gas could result 
in gamma-rays via 7r°-decay. Given that galaxy clusters 
are thought to be capable of accelerating protons to en- 
ergies E > 10 18 eV [193], 7-ray emission could also be 
related to IC emission by secondary electrons generated 
via Bethe Heitler pair production with CMB photons 
(p + 7 — > e + e~ + p) [194,195]. Finally, self-annihilation 
of a dark matter particle (e.g., WIMP) could possibly 
also lead to gamma-ray production [196]. However, de- 
spite these expectations, no extended gamma-ray signal 
has been found so far in VHE observations of local clus- 
ters. The Perseus Cluster, for example, one of the closest 
(d ~ 77 Mpc) and largest clusters has been regarded as a 
particularly promising candidate. This cluster is charac- 
terised by intense X-ray emission and hosts a luminous 
radio halo (size ~ 200 kpc), indicative of the presence 
of a high density of CR protons. Large intra-cluster 
fields (~ 1 — 10/iG) suggest it to be an excellent place 
for particle acceleration. Yet, no (extended) gamma- 



ray excess (but two RG, see above) has been found 
in recent VHE observations (by MAGIC for a total of 
99h in 10/2009-02/2011) [197]. Comparing this with 
simple CR model prediction suggests that the average 
ratio of CR-to-thermal pressure has to be below sev- 
eral percent. Similar (yet somewhat less stringent) CR 
constraints have also been inferred based on the VHE 
non-detection of the Coma Cluster (d ~ 100 Mpc) by 
VERITAS and H.E.S.S. [198,199] (for constraints on CR 
mixing conditions in the Hydra Cluster, see [200]), while 
the VHE non-detection of the Fornax Cluster (d ~ 19 
Mpc) has been used to put constraints on the dark mat- 
ter (velocity- weighted < <jv >) annihilation cross- section 
[201] (cf. also Ref. [300] for recent < <jv > upper limits 
based on central Galactic Halo data). 

Nearby quiescent or "Passive" supermassive Black Holes 
(PBHs), i.e., very massive black holes hosted in galac- 
tic nuclei without bright signatures of broad-band emis- 
sion and with very low luminosity at lower frequency, 
could potentially also produce ultra-high energy cosmic 
rays (UHECR) and VHE gamma-rays by a number of 
magnetospheric processes (see [202] for review). A re- 
cent analysis [203] based on the non-detection of one of 
the prime candidate sources, NGC 1399 (distance^ 20 
Mpc; black hole mass ~ 5 x 1O 8 M ), with current HE & 
VHE instruments allows to put constraints on, e.g., gap- 
type particle acceleration scenarios, suggesting that this 
source is unlikely to be capable of accelerating protons 
to energies beyond 4 x 10 18 eV. 



3 Physics Impact of Recent Results 

3.1 CR and Galactic Gamma- Ray Sources 

Cosmic-rays (CRs) are highly energetic protons, nu- 
clei and electrons which are measured from ~ 10 9 eV 
up to ~10 21 eV with an almost featureless continuous 
spectrum. The most prominent deviation from a pure 
power-law (with index of about —2.8) happens in the 
"knee", where the index changes to ~ —3.1 at about 

4 PeV [204], and the "ankle", a second steepening at 
about 5xl0 17 eV. It is widely believed that CRs be- 
low the knee are accelerated inside our Galaxy. Those 
CRs carry on average as much energy per unit volume as 
the energy density of star light, the interstellar magnetic 
fields, or the kinetic energy density of the interstellar 
gas. One hundred years after their discovery [205], the 
question of their origin is still open. The acceleration, 
accumulation and effective mixture of non-thermal par- 
ticles, through their diffusion and convection in galactic 
magnetic fields, produce the so-called ".sea" of Galactic 
CRs. SNRs are thought to be the main contributor to 
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this CR sea. The main (phenomenological) argument in 
favour of this hypothesis is the CR production rate in 
the Galaxy, W C r ~ (0.3 - 1) x 10 41 erg/s. This rate can 
be easily explained if one assumes that ~10% of the ki- 
netic energy released in supernova explosions (which can 
reach up to 2 x 10 39 erg s _1 kpc -2 ) is devoted to the ac- 
celeration of CRs (see, e.g., Refs. [206-208]). Moreover, 
this highly efficient conversion of kinetic energy of bulk 
motion to relativistic particles can be achieved through 
diffusive shock acceleration (DSA) [7,8,184], which takes 
place in the shocks created in SNRs. A straightforward 
test of the acceleration of CRs in SNRs would be the 
detection of hadronic gamma-rays (TeV 7-rays result- 
ing from hadronic interactions through the production 
and decay of 7r°-mesons) directly from young remnants 
[209] and/or from dense clouds overtaken by the expand- 
ing shells [210]. The DSA theory predicts observational 
features that can be tested with Cherenkov telescopes: 
The characteristic feature of shock acceleration in the 
non-linear regime is a concave shape of the particle en- 
ergy distribution. At low (GeV) energies, the particle 
spectrum is relatively steep with a differential spectral 
index larger than 2, whereas it becomes very hard at the 
highest energies. In the presence of a strongly modified 
shock, the proton spectrum, just before the high-energy 
cut-off, can be as hard as E~ l - b (see, e.g. [24] ). These 
features should then be reflected in the spectrum of sec- 
ondary gamma-rays [28,40,211,212]. However, energy- 
dependent propagation effects of the particles could in- 
troduce significant modifications in the proton spectra, 
in particular in dense regions where most of the gamma- 
rays are produced. This applies to massive molecular 
clouds located outside mid-age SNRs [213,214], as well 
as possible, dense compact condensations inside young 
SNR shells [28]. Such propagation effects can lead to a 
substantial deviation of the observed gamma-ray spectra 
from the parent protons spectrum. 

The new sensitive, wide-field-of-view TeV instruments, 
such as H.E.S.S., have been able to image SNRs in TeV 
gamma-rays, probing the SNR-shell acceleration of ei- 
ther electrons or hadrons up to at least 100 TeV (in case 
of leptonic emission) or a few hundred TeV (for hadronic 
acceleration). The acceleration sites are spatially co- 
incident with the sites of non-thermal X-ray emission, 
strengthening the hypothesis that primary Galactic CRs 
up to the "knee" of the energy distribution are acceler- 
ated in SNRs. However, even if radio and X-ray data sug- 
gest that SNRs are indeed the sources of CR electrons, 
no compelling evidence for the acceleration of protons 
in SNRs has been found up to now, and it is not clear 
whether proton and electron accelerators are of different 
nature. In fact, only a relatively small number of SNRs 
has been detected at VHE (see Sec. 2.1). A possible 
explanation for this is related to the evolution of SNRs, 
which may only be able to accelerate CRs to PeV ener- 



gies in the first ~ 1000 yr, while later the high-energy 
hadrons escape from the system. For instance, a spectral 
cut-off is observed in the bright SNR RX J1713. 7-3946 
at ~4 TeV which implies an exponential cut-off in the 
proton spectrum at ~ 25 TeV. 
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Fig. 20 Radial profiles of 1 TeV gamma-rays calculated for the 
hadronic and electronic scenarios of Ref. [28] in a uniform medium 
(solid) and for a leptonic scenario with unmodified forward shock 
(dashed). The profiles, smoothed with a Gaussian point spread 
function with a = 0.05°, are also shown (thin lines). The trian- 
gles correspond to the azimuthally-averaged TeV gamma-ray radial 
profile as observed by H.E.S.S. From Ref. [28]. 

Still, the production of gamma-rays via pp-interactions 
in dense gas condensations embedded in low-density 
shells represents an interesting scenario which maintains 
a hadronic origin as a viable option with several attrac- 
tive features. The increase in photon statistic in future 
observations with Fermi-LAT and H.E.S.S. should help, 
but may hardly be sufficient to distinguish unambigu- 
ously the contributions of leptonic and hadronic interac- 
tions to the different bands of the gamma-ray spectrum. 
In this regard, CTA has an important role to play. This 
applies, first of all, to a precise measurements of the en- 
ergy spectrum below 1 TeV (down to tens of GeV), and 
above 10 TeV (up to 100 TeV). Morphological studies 
will provide independent and complementary informa- 
tion about the radiation mechanism. The low magnetic 
field, which is a key element in any IC model, allows 
multi-TeV electrons to propagate to large distances, and 
thus to fill a large volume. Because of the homogeneous 
distribution of target photon fields, the spatial distribu- 
tion of the resulting IC gamma-rays will appear quite 
broad. Hadronic models, on the other hand, predict nar- 
rower and sharper spatial distributions, mainly due to 
the enhanced emission in the compressed region of the 
shock, as can be seen in Fig. 20. Using currents instru- 
ments with a typical point spread function of ~ 3 arc 
minutes, this effect is smeared out. An angular resolu- 
tion of 1-2 arc minutes, however, as reachable with CTA, 
would be sufficient to investigate the radial profile with 
enough precision to test these scenarios. 
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A second univocal proof of the origin of CRs in SNRs 
would be the detection of 7-ray emission at extremely 
high energies, the so-called PeVatron (1 PeV=10 15 eV). 
In SNR shocks with relatively low acceleration rates, syn- 
chrotron losses typically prevent the acceleration of elec- 
trons to energies beyond 100 TeV. In addition, at such 
high energies IC emission is suppressed due to Klein- 
Nishina effects. Therefore, a contribution of IC gamma- 
rays to the radiation above 10 TeV is expected to gradu- 
ally fade out. Thus, detecting gamma-rays up to 100 TeV 
would unambiguously establish a hadronic origin of the 
radiation. 

Figure 21 shows an example for the expected X-ray and 
gamma-ray emission of a T = 1000 yr-old proton Pe- 
Vatron calculated for three different distributions of the 
accelerated protons. Both radiation components are ini- 
tiated by the interactions of accelerated protons with am- 
bient gas (of density n = 1 cm -3 ) and magnetic fields 
B = 300 jiG. The gamma-ray emission arises directly 
from the decays of 7r°-mesons, while the X-rays is due to 
the synchrotron radiation from secondary electrons (the 
products of ^-decays). The lifetime of electrons pro- 
ducing X-rays, t syn ~ 1.5 Bj/ 2 (E X /1 keV)" 1 / 2 yr, is 
very short (^ 50 yr) compared to the age of the source. 
Roughly the same fraction of the energy of the parent 
protons is shared between the secondary electrons and 
the gamma-rays from 7r°. However, since the energy of 
the sub- TeV electrons is not radiated away effectively, 
the direct (7r°-decay) gamma-ray luminosity exceeds the 
synchrotron luminosity. The ratio Ly^/L 1 depends on 
the proton spectrum and the injection history, and typi- 
cally does not exceed 0.2-0.3. The noted calculations are 
based on a constant injection rate L p = 10 39 erg/s, giving 
a total energy in protons of W p — L p ■ T ~ 3 x 10 49 erg. 
In order to estimate the X-ray and gamma-ray en- 
ergy fluxes (in units of erg/cm 2 s) from an arbitrary 
PeVatron, the luminosities shown in Fig. 21 need 
to be multiplied with a factor n « 10 _44 (nW p /3 x 
10 49 erg/cm 3 )(d/l kpc) -2 . This indicates that with the 
future CTA instrument, with an anticipated sensitivity 
at 10 TeV of about 10 -13 erg/cm 2 s (cf. Fig. 1), all 
Galactic sources capable of accelerating particles to PeV 
energies with nW p ^ 10 49 erg/cm 3 should become de- 
tectable up to distances of 10 kpc. 
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Fig. 21 Broad-band radiation of a PeVatron initiated by inter- 
actions of protons with the ambient gas. The curves are based 
on three different proton distributions: a single power-law with an 
exponential cutoff (E~ a exp(-E/ Eq)P , with a = 2, E = 3 PeV) 
but different indices (3 = 1 (solid curve) and /3 = 1/3 (dashed 
curve), respectively, and a broken power-law where the spectral 
index changes at E = 1 PeV from a = 2 to a = 3. From Ref. [23]. 

With respect to PeVatron considerations, the relevance 
of CR escape from the shell confinement has become 
more evident in the last years. Efficient acceleration 
and confinement of mufti- TeV particles in SNRs can last 
less than several hundred years after the explosion, see 
e.g. [215], effectively constraining the number of SNRs 
emitting as PeVatrons. However, these escaping CRs 
could still be observable if certain conditions in the sur- 
rounding environment are fulfilled, such as the presence 
of massive molecular clouds which could provide a dense 
target to trap those running-away CRs for hadronic in- 
teractions, before they diffuse away and get integrated 
in the CR sea. Depending on the distance between the 
SNR and the molecular cloud (of typical mass of 1O 4 M ), 
the time of particle injection into the interstellar medium 
and parameters such as the diffusion coefficient, different 
gamma-ray spectral shapes might be expected, reach- 
ing TeV flux levels detectable by current instruments 
[213,216]. 

Finally, the search for CRs accelerators also needs to be 
extended to other source classes. Different types of accel- 
erators have been proposed to contribute to the Galactic 
CRs, such as superbubbles, see e.g. [217-219], or the rem- 
nants of gamma-ray bursts in our Galaxy [220]. More- 
over, the CR flux measured at Earth might deviate from 
the Galactic one, if it would be dominated by a single 
nearby source (or a few local sources) in the surrounding 
of the solar system. The recent CR "anomalies" dis- 
covered by the PAMELA experiment, such as the very 
high content of positrons in their leptonic component 
[221], or the significant differences in the energy spec- 
tra of protons and alpha-particles [222], might be due to 
an oversimplification of a CR single-class-of-accelerator 
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origin. In fact, a recent analysis [223] has reported vari- 
ations of the CR spectrum in nearby molecular clouds 
embedded in the Gould Belt, which may indicate such 
deviation. Constraints on the CR diffusion coefficient in 
clouds nearby to SNRs also indicate a lower value than 
the one measured from the composition of the local CR 
spectrum (see, e.g. [224]), suggesting that the CR diffu- 
sion coefficient might not be isotropic, which could lead 
to differences in the CR spectrum in different regions of 
the Galaxy. 

In this concern, the Galactic plane survey planned with 
the CTA observatory could result in the discovery of new 
and exciting classes of PeVatron gamma-ray sources of 
unknown origin [22]. 

3.2 Relativistic Outflows and AGNs 

Powerful relativistic outflows and jets, with speeds ap- 
proaching that of light, have been seen in a number of 
astrophysical sources including Active Galactic Nuclei 
(AGNs), Pulsars and Microquasars (MQs). If these out- 
flows point towards us, as thought to be the case in e.g. 
blazars, relativistic effects can dramatically change their 
appearance. The dynamics of and the non-thermal pro- 
cesses in these outflows are important fields of active 
research. Some examples of current frontiers relevant to 
TeV Astronomy are mentioned below: 

3.2.1 Relativistic pulsar winds: 

Outside the light cylinder, pulsars are believed to launch 
highly magnetised, electron-positron winds with termi- 
nal flow speeds possibly reaching Lorentz factors up to 
T w ~ 10 6 (for recent review see e.g., [225-227]). At the 
location where the interaction with the ambient medium 
(SNR or ISM) causes the wind to slow down, a "termina- 
tion" shock is formed. It is here where efficient particle 
acceleration is thought to energise the particles radiat- 
ing in what is observed as pulsar wind nebula (PWN). 
The theory of pulsar winds is usually based on the MHD 
description of relativistic outflows, with its electromag- 
netic structure being strongly dependent on the current 
it carries. In the simplest (force- free) approximation the 
plasma flow in the relativistic case becomes asymptoti- 
cally radial. The magnetic field is considered frozen into 
the flow, so that the field structure in the far wind zone 
becomes almost purely azimuthal/toroidal. The pulsar 
wind is commonly believed to be "cold" (i.e., charac- 
terised by ordered motion only; in the co-moving frame 
all particle are stationary) and therefore mostly unob- 
servable (e.g., there is no synchrotron radiation). Bulk 
inverse Compton scattering of ambient photons, how- 
ever, could still reveal its presence ([228,229]; see also 
below). In the current MHD picture (cf. Fig. 22), the 
wind is launched as a highly magnetically (Poynting flux) 
dominated flow with a = B 2 / A-Knm e c 2 Y w ^> 1 (ratio of 



Poynting to kinetic energy flux in the unshocked wind). 
Spectral modeling for the Crab Nebula suggests that 
close to its termination shock most of this energy should 
have been converted to the kinetic energy of the bulk mo- 
tion (i.e., acceleration of the wind), resulting in a <C 1 
[230]. How and where exactly this is achieved, however, 
is not yet completely understood ("a-problem"). In an 
unconfined (radial) ideal MHD outflow the asymptotic 
Lorentz would eventually become comparable to that of 
the fast magnetosonic mode T WiOQ ~ ^Jo only, implying 
~ (a r ) 2 / 3 > 1 (e.g., typically a ~ 10 4 , T ~ 100), 
i.e. the wind would stay Poynting-dominated. If this 
were to be true, efficient first-order Fermi acceleration 
at the termination shock would not be possible (due to 
a low compression ratio and a relativistic downstream 
speed). 

The recent detection of pulsed gamma-ray emission be- 
tween - (100 - 400) GeV from the Crab pulsar [57,58] 
now provides for the first time strong observational sup- 
port for the presence of a (cold) ultra-relativistic pulsar 
wind and indicates that wind acceleration to T w > 10 5 
could already occur very close to the pulsar, i.e., at 
distances between (20-50) light-cylinder radii) [60], see 
also Fig. 4. Inverse Compton up-scattering (in the KN 
regime) of the pulsed magnetospheric X-ray emission by 
pulsar wind electrons seems to be responsible for VHE 
emission up to ~ T w m e c 2 = 500 (T w /5 x 10 5 ) GeV. 




Fig. 22 Anticipated pulsar wind structure in the case of the 
Crab. From Ref. [60]. 

The apparent rapid conversion of Poynting to kinetic en- 
ergy flux is certainly not easily accommodated in most 
of the current scenarios. It could be that this conver- 
sion attests to a longitudinal current distribution that 
is smaller than the Goldreich- Julian one [231]. On the 
other hand, if the pulsar wind would still possess a 
non- negligible magnetisation (a ^ 0.1) approaching the 
termination shock, cross-field transport would be sup- 
pressed, making efficient operation of first-order Fermi 
particle acceleration unlikely. In such a case, excitation 
of fast growing micro-turbulence (with 5B/B > 1 on 
scales X < tl) would be required, but seems rarely to oc- 
cur under the noted conditions [232]. Magnetic field an- 
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nihilation by shock-driven reconnection in a striped wind 
[233-235] could perhaps offer an escape route, both by 
producing "non-thermal" (power-law type) particle dis- 
tributions in the reconnection electric field, and by gen- 
erating sufficiently strong turbulence for further Fermi 
acceleration. Alternatively, other non-ideal MHD effects 
(such as conversion of the wind from a sub-luminal MHD 
into a superluminal electromagnetic wave, expected to be 
accompanied by a substantial energy transfer from fields 
to particles [236]) may become important in outer wind 
regions. Similar challenges (perhaps even more severe 
given the fact the here the termination shock is located 
much closer to the pulsar) are encountered in the gamma- 
ray binary context (cf. 3.2.3). Further VHE studies like 
the ones noted above may soon help us to settle a fun- 
damental issue in current pulsar wind theory. 

3.2.2 Relativistic jets in AGNs: 

AGNs are observed to exhibit collimated (often two- 
sided) relativistic outflows extending from sub- up to 
hundreds of kilo-parsecs. Radio VLBI/VLBA observa- 
tions of blazar jets have revealed significant apparent su- 
perluminal motion (f3 a = v a /c > 1) of individual jet 
components on parsec-scale propagating away from the 
core. The velocity distribution of the fastest measured 
radio jet components peaks at f3 a ~ 10, but seems to 
possess a tail extending up to ~ 50 [237,238]. This in- 
dicates that characteristic flow speeds in AGN jets may 
reach bulk flow Lorentz factors up to T ~ @ a /fi ~ 50. 
Current evidence, however, suggests that the apparent 
(parsec-scale) radio jet speeds in the TeV-detected HBL 
sources are much slower than the speeds in the above 
noted, radio-selected blazars [239]. In fact, the appar- 
ent lack of superluminal speeds in the parsec-scale ra- 
dio jets of TeV-HBLs seems to suggest that the bulk 
Lorentz factor in these object is only modest (r ~ 2 — 3), 
in contrast to estimates of the (sub-parsec-scale) bulk 
Lorentz factors based on simple SED modelling. This 
could attest to the presence of some velocity gradient 
in the flow, such as a (longitudinally) decelerating flow 
[241,242] or a (transversally) structured jet with a fast 
moving spine and a slower moving sheath [242]. If this 
were to be true, the radiative interplay between differ- 
ent zones and the possibility of further particle acceler- 
ation [243], could easily lead to complex emission pat- 
terns, thereby severely limiting classical inferences based 
on one-zone synchrotron- Compton modelling. Note that 
a spine-sheath (fast flow/slow flow)-type topology may 
be naturally expected within the context of current MHD 
models for the formation of relativistic jets. In principle, 
these models rely for power extraction on a source of ro- 
tation, i.e., either on the accretion disk (in a Blandford- 
Payne scenario) or the spin of the black hole itself (in 
a Blandford-Znajek scenario), see e.g. [244,245] for a 
recent illustration. The extent to which one of these sce- 



narios dominates in AGNs is not fully understood [246]. 
There are recent indications that the black hole spin is 
less important for powering AGN jets than previously 
believed [247], and that the radio power output in FR I 
type sources, believed to be the parent population of BL 
Lacs, is consistent with a rather modest spin distribution 
[248]. Perhaps most likely a mixture occurs in which a 
pair-dominated, fast-moving (r ^ 10) ergospheric-driven 
jet is encompassed by a slower- moving (r £ 10), disk- 
driven electron-proton jet. Clearly, proper modelling of 
the VHE source states in TeV blazars can shed light on 
these issues by constraining the required distribution of 
bulk flow Lorentz factors and the required power content 
in non-thermal particles. 

The observed variability in VHE blazars with evidence 
for sub-hour variability (cf. Table 2) in several sources 
points to a compact VHE emitting region of size R ^ 

L>cA£ V he(1 + z) ~ 10 15 (L>/10)(A£ V HE/lh) cm ' As ~ 
suming a fiducial jet opening angle i ~ 1/r, as often 

done, this would (without further considerations) trans- 
late into distances of d ~ YR <C 1 pc. While this often 
appears to be a good working assumption, such an infer- 
ence seems severely challenged in the case of the distant 
VHE blazar PKS 1222+216 (4C 21.35, z = 0.432, BH 
mass of a few times 1O 8 M ). PKS 1222+216 is a bright 
FSRQ where the strong disk UV continuum (Ld ~ 10 46 
erg/s) is thought to photo- ionize what is called to be 
the "broad line region" (=BLR; often taken as a quasi- 
spherical ensemble of gas clouds moving with speeds of 
several 1000 km/s at distances of d ~ 0.1 pc). If the VHE 
gamma-rays were to originate inside the BLR, strong 77- 
absorption features (due to pair production in the dense 
BLR UV photon field) should become apparent in the 
observed VHE spectrum. Despite this expectation, how- 
ever, the (EBL-deabsorbed) VHE spectrum as seen by 
MAGIC in June 2010 (coinciding with a high GeV source 
state) seems well compatible with a single power-law of 
photon index 2.7 ± 0.3 between 3 GeV and 400 GeV, 
without evidence for an intrinsic cut-off [249], see Fig. 
23. The obvious inference that the VHE emission in this 
source may simply originate from beyond the BLR has 
spurred considerable theoretical efforts and led to a rapid 
proliferation of models, e.g. [251-254]. One possibility 
is that efficient gamma-ray production occurs when the 
parsec-scale jet experiences a strong re-collimation shock 
[250,251]. 
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Fig. 23 Differential energy spectrum of FSRQ source PKS 
1222+216 as measured by MAGIC. Upper limits are shown as black 
arrows. The solid line is the best fit of the measured spectrum to a 
power law. From Ref. [249], reproduced by permission of the A AS. 



Alternatively, the VHE gamma-rays could also result 
from energetic, secondary leptons produced in photo- 
hadronic (n+7 — >> ir) and 77-interactions of neutrons and 
UHE 7-rays with infrared photons of the more distant 
dust torus, provided efficient UHE proton acceleration 
in the inner jet, generating an escaping beam of UHE 
neutrons via photo-hadronic processes, and Doppler fac- 
tor D ~ 100 do occur [252]. Other models invoke oscil- 
lations of VHE source photons into axion-like particles 
[293] or magnetic reconnection events [254] . What seems 
to be in need of further study at this stage, however, 
is the extent to which more a realistic BLR modelling 
could perhaps allow to relax some of the opacity con- 
straints (see [255] for a first step). As the BLR is not 
spatially resolved on direct images, reliable information 
on its detailed structure and kinematics is usually not 
straightforwardly available. It seems likely, for example, 
that in reality the BLR is more complicated, e.g., [256- 
257], being composed of multiple (geometrically differ- 
ent) regions (e.g., an inner, disk- ionized and an outer, jet- 
energized part; spherical or axial symmetric etc) and that 
the inferred BLR velocity is a superposition of different 
components (turbulence, in- and outflows, rotation etc). 
Interestingly, the H^ emission line flux in PKS 1222+216 
seemed to have varied rather little during the past years 
(within a factor ^ 1.7), despite relatively large (factor 
~ 4.3) optical continuum variations [259,260], which may 
indicate that the overall BLR is not strongly affected by 
the observed variable and supposedly highly beamed jet 
emission. Optical and VHE studies along this line could 
obviously help to shed more light on the BLR structure 
in powerful bazars. 

For M87 (being more massive than PKS 1222+216 by 
an order of magnitude) the observed day-scale VHE 
activity is already comparable to the light travel time 
across Schwarzschild-radius distances. Under the condi- 



tion of moderate Doppler boosting (D < 3, as expected 
for a misaligned jet-source) the observed variable and 
hard VHE emission should originate from near to the 
black hole. This is facilitated by the fact that M87 
is highly under-luminous (compared to Eddington, and 
in contrast to PKS 1222+216), making it possible for 
VHE 7-rays to escape from its inner jet environment. A 
multitude of VHE emission scenarios, relying on either 
hadronic and/or leptonic radiation processes, have been 
developed and applied within recent years to get insight 
into the real nature of the source (see [170] for review 
and references). In the context of leptonic models, one 
interesting variant assumes relativistic (Petschek-type) 
reconnection as acceleration mechanism of electrons to 
high energy and external Compton scattering to account 
for the VHE emission [261,262]. Efficient reconnection 
could possibly lead to an additional relativistic veloc- 
ity component of the ejected plasma (in the bulk jet 
frame), making differential Doppler boosting possible, 
thereby circumventing the problem of a modest Doppler 
factor for the general bulk flow. It is, however, not clear 
at the moment whether the required conditions (e.g., 
high magnetisation for an electron-proton jet, negligi- 
ble guide field) are likely to prevail in AGN jet. Within 
hadronic models, on the other hand, one recent develop- 
ment utilises VHE gamma-ray production via inelastic 
proton-proton (pp) collisions, assuming that interactions 
of a red giant star or a massive, dense gas cloud with the 
base of the jet could occasionally introduce the required 
high amount of target matter [263]. Recent modelling 
suggests that such a scenario could possibly account for 
the observed temporal and spectral properties if the jet 
is powerful enough and a sufficiently large fraction of 
it could be channeled into VHE gamma-ray production. 
An alternative explanation has been pursued in the con- 
text of magnetospheric models (see [202] for review and 
references), motivated by the observed day-scale VHE 
variability (comparable to Schwarzschild scales) and the 
observed radio- VHE connection in 2008 [175]. In a re- 
cent scenario [264], pair creation in a hot accretion flow 
(ADAF) is considered to lead to the injection of pri- 
mary electrons and positrons into magnetospheric vac- 
uum gaps, in which they are quickly accelerated to very 
high Lorentz factors. Curvature emission and, in par- 
ticular, Compton up-scattering of (ADAF) disk photons 
then results in 7-rays with a spectrum extending up to 
10 4 TeV. For low accretion rates, photons with energies 
below several TeV could escape 77-absorption. On the 
other hand, IC photons with energies above ~ 10 TeV 
could interact with the ambient radiation field initiating 
pair cascades just above the gap. This would lead to a 
large pair multiplicity and ensures that a force-free out- 
flow becomes established, accounting for the jet feature 
as seen in radio VLBA. 

Obviously, due to its proximity and low luminosity, M87 
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seems to bear the rare potential to give us important 
insights into the near black hole environment of active 
galaxies. Future, more sensitive VHE observations will 
be important to study the spectral variability and fur- 
ther constrain the time scales of the VHE flux variations, 
while simultaneous radio observations will have the po- 
tential to pin down the location of the VHE emission on 
spatial scales ^ 10 2 r s . 

Perhaps the most prominent example for extreme short- 
time VHE variability is the BL Lac object PKS 2155- 
304. In July 2006, clear evidence for minute-timescale 
(t v ~ 200 sec) VHE variability was found during a dra- 
matic outburst (with flux level varying between 1 and 
15 Crab units) [137]. The remarkable data set made a 
detailed temporal characterisation of the flaring period 
possible, showing that the variability power is distributed 
over (temporal) frequencies according to a red- noise-type 
PSD (power spectra density) oc with f3 ~ 2 in the 
frequency range (10~ 4 - 10~ 2 )Hz, see Fig. 24 [137,265] 
(cf. also [266] for a related PSD GeV variability analysis 
below - 10~ 6 Hz). 
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Fig. 24 Power spectrum density (PSD) of the light curve of PKS 
2155-304. The grey shaded area corresponds to the 90% confidence 
interval for a light curve with a PSD oc v~P . The dotted horizontal 
line shows the average noise level. From Ref. [137]. 

Similarly to (previous) results in the X-ray domain, a 
correlation between (absolute) rms (root mean square) 
variability amplitude and mean VHE flux has been 
found. This is known to be characteristic of a non-linear, 
log-normal stochastic process and suggests that the un- 
derlying process driving the VHE variability is of mul- 
tiplicative (not additive) nature (see [267], but cf. [268] 
for a possible caveat). The extreme short-term variabil- 
ity has attracted considerable attention, with a number 
of interesting models developed to account for it (e.g., 
spine-sheath [269], reconnection/ mini-jets [270], jet-star 
interactions [271], non-ideal MHD effects [272] or tur- 
bulence [273]). Like in Mkn 501, the VHE variability in 
PKS 2155-304 seems to occur on timescales much shorter 



than the horizon light-crossing time for a typical black 
hole mass inferred from its host galaxy luminosity. This 
could suggest that the variability originates in a small 
region within a highly relativistic outflow having bulk 
Lorentz factors ~ 50 [274] (but see [275] for an inter- 
esting alternative based on time-dependent modelling). 
The detailed temporal characteristics of the flare (PSD 
and log-normality) have not yet received as much atten- 
tion. This situation needs to be remedied as it could offer 
the potential to distinguish between the different models 
proposed. Linear RMS-flux relations, for example, are 
known from black hole X-ray binaries and believed to be 
an essential property of the accretion flow around black 
holes [276], likely produced by small, independent fluc- 
tuations in the accretion rate [277]. An application of 
such a scenario to PKS 2155-204 is possible, but seems 
to require the presence of a supermassive binary black 
hole system [278]. The increased sensitivity of future in- 
struments like CTA will allow us to probe the temporal 
characteristics of flaring TeV blazars beyond the usual 
minimum variability considerations and thereby offer an 
important diagnostic tool to uncover the real nature of 
the source. 

3.2.3 Relativistic outflows in galactic VHE binaries: 
Galactic VHE binaries are composed of a compact ob- 
ject (a stellar mass black hole in accreting binary sys- 
tems=Microquasar [MQ] scenario as in e.g., Cygnus X- 
1; or a ms-pulsar in pulsar wind binaries as in PSR 
B 1259-63) and a bright companion star orbiting around 
each other. Complex interactions between the stellar 
wind and the outflow produced by the compact object 
(a mildly relativistic MQ jet, or an ultra-relativistic pul- 
sar wind) are naturally expected to occur, potentially 
triggering efficient dissipation and gamma-ray produc- 
tion through a variety of processes (see [279] for review). 
Current research frontiers in this context include the fol- 
lowing: 

(1) The interactions of a relativistic pulsar wind with the 
wind from a stellar companion can result in the formation 
of an extended cometary-shaped radio structure rotating 
with the orbital period [280]. Electrons are thought to 
be efficiently energised in the termination shock region, 
producing HE and VHE gamma-rays by inverse Comp- 
ton up-scattering of the stellar UV photons whilst being 
advected away [280,281]. In the case of the pulsar wind 
binary PSR B1259-63, observational evidence for vari- 
able and extended (AU-scale) radio emission has been 
found, demonstrating that extended (jet-like) structures 
can be produced in pulsar wind binaries [282]. Rotating 
jet-like radio features have also been detected in LS 5039 
and LS I +61 303 [283-285], whose natures are not yet 
established. In analogy to PSR B 1259-63, these findings 
could suggest a pulsar wind binary interpretation. On 
the other hand, jet-stellar outflow interactions may be 
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capable of producing similar features [286,288], and MQ 
jet precession (with variable Doppler boosting) has been 
recently favoured for LS 5039 based on the finding that 
its radio emission displays double-sided structures chang- 
ing to one-sided ones [284]. The situation could be even 
more complicated in that the extended radio emission 
may partly also be related to synchrotron emission from 
secondary pairs generated in 77-interactions [287]. What 
seems to be an open issue for a pulsar wind interpretation 
when applied to LS I +61 303 is that for a reasonable 
7-ray production efficiency the required PW momentum 
flux significantly exceeds the stellar (Be) wind one, mak- 
ing the formation of elongated cometary-shaped radio 
structures unlikely [288]. 

(2) The location of the VHE emitting region is not yet 
understood and different (not mutually exclusive) emit- 
ting sites are conceivable (see also above). In principle, 
VHE gamma-rays could be produced inside the binary 
(pulsar wind zone or jet base), on binary scales (wind 
collision/interaction region) or on even larger scales (jet 
termination). In the case of LS 5039, the VHE data 
favour an origin of the TeV emission on binary (~ 10 12 
cm) scales [90,289]. The possibility of efficient particle 
acceleration (at shocks) then becomes an open issue. If 
particle acceleration occurs inside a MQ jet, only mildly 
relativistic flow speeds may be expected. If, on the other 
hand, particle acceleration is believed to occur in the 
shock region formed by the encounter of the pulsar wind 
and the stellar wind (standoff distance much closer to the 
pulsar compared to the termination shock for isolated 
pulsars), efficient first-order Fermi acceleration may be 
difficult to achieve (cf. section 3.2.1.). A comparison of 
the HE and VHE emission in PSR B1259-63 indicates 
that the spectacular GeV flare seen by Fermi is not ac- 
companied by a similar flare in the TeV regime, suggest- 
ing that the GeV and the TeV emission are produced 
by different mechanisms [290]. Similarly, in the case of 
LS 5039 and LS I +61 303, the exponential cut-off at a 
few GeV, combined with the hard spectrum at TeV and 
the anti-correlation between the two bands, points to- 
wards different origins for their GeV and TeV emissions 
[291]. Detailed numerical simulations of stellar/pulsar 
wind or stellar wind/MQ jet interactions on different 
scales, see e.g. [292-294], are expected to improve our 
understanding of the dominant VHE emitting region. 

(3) Finally, interactions of the jet with the ambient 
medium/ISM (termination) may also result in efficient 
particle acceleration facilitated by, e.g., re-collimation, 
forward or reverse shocks, and possibly lead to VHE 
gamma-ray production that may eventually become de- 
tectable [295]. If so, then this could allow to more tightly 
constrain the jet kinetic power and the magnetic field in 
the interaction region. 

While further high-resolution radio observations bear the 
promise to reveal the true origin of the changing radio 



morphology in gamma-ray binaries, high(er) sensitivity 
VHE observations with dense temporal coverage are ex- 
pected to shed new light on the non-thermal mecha- 
nisms responsible for TeV gamma-ray production, and 
will eventually help to disclose the real nature of these 
sources. 



4 Conclusion and Perspectives 

The remarkable achievements of observational gamma- 
ray astronomy over the last decade, and the recent theo- 
retical and phenomenological studies of acceleration and 
radiation processes in astrophysical environments, fully 
justify the further exploration of the sky at high and 
very-high energies. Although generally the main moti- 
vations of gamma-ray astronomy remain unchanged, the 
recent observational results have revealed new features 
which in many cases require revisions of existing theo- 
retical models and formulations of new concepts. It is 
expected that over the next decade the ongoing opera- 
tion of Fermi will be accompanied by observations with 
the current (H.E.S.S., MAGIC, VERITAS) and planned 
(CTA, HAWC) ground-based detectors. The data ob- 
tained in the enormous energy range from 100 MeV to 
1 PeV will provide very deep insight into a number of 
problems of high-energy astrophysics and fundamental 
physics. 
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